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ABSTRACT 


This report covers the scope of Mod. 4 of Contract NAS8-28358, covering 
a period of less than one year which is interior to the total time span of 
the contract. Thus what follows is a "final report" only in the sense of 
providing final documentation of those tasks which have been completed during 
this period. Specifically, this report includes the following subjects: 

"The Influence of Spacecraft Flexibility on System Controllability and 
Observability;" "Commutativity of Coordinate Truncation and Transformation 
Matrix Inversion for Flexible Spacecraft Dynamic Analysis;" and "Matched 
Asymptotic Expansion Modal Analysis of Rotating Beams." 
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OVERVIEW 


Contract NAS8-28358 began on 15 February , 1972, and is currently 
scheduled through 14 December, 1974. The present report covers only' "Mod. 

4" of the contract, which was awarded in late Spring of 1973 and covers 
the period from 15 February 1973 to 15 February 1974. In this report there 
appears a final documentation of those tasks presently deemed completed; 
other work still in progress will be documented at the conclusion of "Mod. 

5" of this contract, in December of 1974. 

The three topics treated here are as follows: 

Chapter 1 . "The Influence of Spacecraft Flexibility on System Con- 
trollability and Observability;" 

Chapter 2 . "Commutativity of Coordinate Truncation and Transformation 
Matrix Inversion for Flexible Spacecraft Dynamic Analysis;" 

Chapter 3 . "Matched Asymptotic Expansion Modal Analysis of Rotating 
Beams." 

These chapters have been written by the Principal Investigator in such 
a way that each can be submitted (perhaps in abridged form) for conference 
presentation and/or journal publication. 

Other topics under study during the Mod. 4 period of the contract 
include : 

(1) Development of general-purpose simulation equations for arbitrary 
spacecraft; 

(2) Evaluation of the influence of sensor and actuator location on 
flexible spacecraft performance; and 

(3) Preliminary evaluation of the concept of the disturbance-insensitive 

control system. 
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During the contract period, support funds were extended to Dr. Yoshiaki 
Ohkami and (to a minimal degree) Mr. Oluyemisi Olusola, a student who is no 
longer with the project; substantial contributions to total progress have 
also been made by Mr. Robert Skelton and Mr. Stewart Hopkins (both current 
students) and by Dr. Andre Colin (a former student). 
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CHAPTER I 


The Influence of Spacecraft Flexibility 

A 

on System Controllability and Observability 

ti 

Y. Ohkami 
and 

P. Likins^ 

ABSTRACT. Literal criteria are developed for the controllability and observa- 
bility of general models of flexible spacecraft. Results are interpreted in 
special cases and in physical terms, permitting in some cases the identifica- 
tion of uncontrollable and unobservable states simply by examination of scalars 
composed of modal parameters and location matrices for sensors and actuators. 

A procedure is established for isolation of uncontrollable states, whereby 
sensor and actuator configurations assure that uncontrollable flexible mode 
states are also unobservable; in many applications such states can then be re- 
moved by coordinate truncation. 


Based on research supported by NASA Marshall Space Flight Center under Contract 
NAS8-28358. 

* 

Professor, School of Engineering and Applied Science, University of California, 
Los Angeles, California, USA. 

NASA International Fellow on leave at the University of California, Los 
Angeles, from the National Aerospace Laboratory, Tokyo, Japan. 
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INTRODUCTION 


Scientific spacecraft now in the planning stages (as typified by the NASA 
Large Space Telescope) are vehicles of dynamically significant flexibility 
which must be capable of maintaining inertial orientations to within thousandths 
of an arcsecond over periods of hours. Sensors and actuators are typically 
distributed over the spacecraft, and influenced in their performance by vehicle 
flexibility. The attitude control problems posed by these spacecraft are of 
unprecedented difficulty, and their solution will require unparalleled preflight 
investigation and perhaps conceptual Innovation, such as system optimization or 
adaptive control. Although problems of spacecraft attitude control have in the 
past been resolved primarily by a combination of digital computer simulation 
and the traditional design procedures of linear control theory, in the future 
the techniques of modern control theory must increasingly be employed. The 
purpose of this paper is to examine the fundamental concepts of controllability 
and observability as they apply to flexible spacecraft, in order to establish a 
foundation on which to build future applications of modern control theory. As 
our results emerge, it will become apparent that they are useful not only in 
projected control system synthesis but also in the modal coordinate truncation 
process which is essential to system simulation. 


i Preceding page blank . 
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EQUATIONS OF STATE 


Although a variety of mathematical models have been adopted for flexible 
spacecraft, the most common idealization consists of a rigid primary body with 
attached elastic appendages, which are themselves modeled sometimes as continua 
[1] but most often as discretized assemblages of particles or rigid nodal bodies 
interconnected by massless elastic elements called finite elements. [2 J In 
some cases, the appendage mass is distributed throughout the finite elements. [3] 
Whether the original appendage model is continuous or discretized, its displace- 
ments relative to the primary body are finally characterized by distributed or 
modal coordinates, whereas the inertial orientation of the primary body is char- 
acterized by discrete coordinates, such as a set of three attitude angles. In 
many applications additional discrete coordinates are employed to describe the 
motion of system components, such as rotors and scientific instruments. A rep- 
resentation of vehicle kinematics in terms of a combination of discrete and dis- 
tributed coordinates is called a hybrid coordinate formulation. 

In what follows, it is assumed at the outset that the vehicle model consists 
of a rigid primary body with discretized elastic appendages having n rigid nodal 
bodies. After coordinate transformation, the resulting dynamic equations have a 
structure that permits their interpretation as representative also of appendage 
models consisting of elastic continua or distributed-mass finite element sys- 
tems. An example in Appendix A Illustrates the applicability of our results to 
continuous appendages. 

It is further assumed in this paper that the vehicle experiences only small 
deviations from a nominal state in which the central body is of inertially fixed 
orientation and the appendages have no deformations. Moreover, the vehicle has 
no inherently nonlinear elements, so that all analysis can be based on linearized 
variational equations. Implications for nonlinear systems can sometimes be 
drawn from the physical interpretations of the conclusions presented here. 
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( 1 ) 


From Eq. (278) of [2], the vehicle rotational equations take the form 

T.. 

10+ Aq - T 

where T is the 3x1 matrix representing the external torque vector about the 

vehicle mass center c f for an orthogonal vector basis in the primary body b, 

* 

1 is the 3x3 inertia matrix of the total vehicle for c, 0 is the 3x1 matrix 
1-2-3 inertial attitude angles of b, q is a 6nxl matrix consisting of a se- 
quential ordering of 6x1 matrices typified for nodal body j by [u^ 3^ 3^ 

where for ot = 1,2,3, u^ is the translation of the nodal body mass center relative 
to b in the direction defined by axis a fixed in b, and 3^ is the small rota- 
tion of nodal body j about axis a. The 6nx3 matrix A is established by the 

geometry and mass distribution characteristics of the appendages, as in Eq. (278) 
of [2]. - 


The appendage deformation equations in Eq. (277) of [2] may be written 
M'q + Kq + A9 = X 

under the assumption that the resultant external force applied to the vehicle is 
zero. Here K is the appendage stiffness matrix and M* is the inertia matrix 
appropriate for the appendage attached to a translationally free but rotation- 
ally constrained rigid body. Both M' and K are symmetric and positive definite, 
and all quantities in Eqs. (1) and (2) are real. The 6nxl matrix X consists 


of n 6x1 partitions for the n nodal bodies, 
Cfi f 2 f 3 *1 ^2 where for a = 1.2,3, 

ex 


as typified for the j th body by 
and are the scalar components in 


direction a of the external force and torque applied to nodal body j . in the 


present application, X contains only forces and/or torques associated with atti 


tude control actuators. Although in practice these actuators might themselves 
be defined by their own dynamical equations, it serves our present objectives to 
treat them as pure thrusters or torquers. 
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The relationship between X and T can be expressed in general terms involv- 


ing the torque applied directly to b and the forces and torques in X as 

T = <X + T o ( 3 ) 

for some 3x6n matrix K which sums the torques on the nodal bodies and cross- 
multiplies the forces by the vector locating nodal body j relative to c; thus 
1 « ‘ 


< = [ft ! u 3 J ft ! u 3 J — ! a n ! u 3 ] 


(4) 


where U 3 is the 3x3 unit matrix and in terms of the scalar components , R-’ s rJ 

12 3 


of R J in the vector basis fixed in b. 


[ft j ] A 





(5) 


In Eq. (3) there are potentially 6n+3 independent controls embodied in T and X. 

o 

Ift the simplest case the 3 elements of T are each identified as independent 
controllers (each perhaps associated with a single actuator, such as an idealized 
control moment gyro); then it is most convenient to replace X in Eq. (2) by 

A = ^c T (6) 


in which the 6nx3 matrix establishes the location and type of the attitude 

control actuators. For example, if the actuator is conceived as a pure torquer 
t h 

located on the k nodal body, then 


K ■ |° | o|— j° J “ 3 ; — | ° i ° ] (7) 

where U 3 is a 3x3 unit matrix appearing in the 2k th position of# . If instead 
the actuation system consists of a system of synchronized attitude control jets 
responding to only three independent commands, then Eq. (6) still provides an 
appropriate actuation model, although#^ becomes more complex than the simple 
case illustrated by Eq. (7). (See the first example, case (c), in Appendix A.) 
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Since the nature of the' three-axis control problem clearly demands at least 
three independent controllers, Eq. (6) represents a minimal actuator model. 

In the next section we address the question: Is the system controllable 

with the minimal actuator model? An affirmative answer implies controllability 
with the more general actuator model of Eq. (3), but the system could theoret- 
ically be uncontrollable with the minimal model of Eq. (6) and yet controllable 
with the model of Eq. (3). 

The dynamical equations (1) and (2) must be augmented by an observation 
equation. Again we are concerned only with the minimal sensor model, which 
provides only three attitude readings. Thus we adopt the observation equation 

y = 0 +& T q q (8) 

where the 6nx3 matrix ^ locates the sensors* For example, if a sensor on the 
k th nodal body provides the three inertial attitude angles of that body, then 
is identical to the matrix illustrated in Eq (7), and the sensor output is 

y = 0 + 0 k (9) 

Physical arguments suggest that observability is unchanged by the acquisi- 
tion of attitude rate information at the same locations for which attitude infor- 
mation is established by Eq. (8). Generalization of this minimal sensor model 
to permit more than three independent attitude measurements may, however, ex- 
tend observability to additional states. 

In the hybrid coordinate formulation, Eqs. (2) provide the basis for the 
transformation 

q “ ^ (10) 

in which <{> is a 6nx6n matrix whose columns are the eigenvectors of that portion 
of Eq. (2) involving only q. The noted properties of M' and K assure that these 
eigenvectors are independent, and that (j) is nonsingular and has orthogonality 
properties such that 4> M'<}> and are diagonal. (See [2], pp, 43 and 53.) 
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Eqs. (1), (2), (7), (8), and (10) become, after pre-pultiplication of Eq. 

(2) by <() T and of Eq. (1) by I \ 

6 - Jri = u (11) 

n + a 2 n - 69 = (<()' r M , (j))“ 1 4) T i^ c i ,k u ( 12 ) 

y = 6 + sfyr) (13) 

where 

*-l 

u A I T 
6 A -(<j> T M'<j>) -1 4> T A 

T T T 

J A -I A cf> = I 6 (<J>‘ V<()) 

2 T -IT 

a = (*V*) 

Eqs. (11) - (13), although obtained here for an appendage model [2] in 
which all mass is concentrated in nodal bodies, apply also to vehicles with con- 
tinuous appendages [1] or with distributed -mass finite element appendage models 

[3] . In any case the dimension of n represents the number of appendage modes, 

which henceforth is called N (rather than the symbol 6n which is appropriate 

only for the fully discretized model without coordinate truncation.) 

T 

In what follows the eigenvectors in 4> are normalized so that <f> M'<fi = U^; 

2 2 

then the diagonal elements o ^ of the NxN matrix o represent squares of non- 
zero natural frequencies of vibrations which could occur independently of each 
other if the primary body were translationally free but constrained against 
rotation, as noted in [2], page 55. 

In terms of the (2N + 6)xl state variable 


Eqs. (11) - (13) may be written as 


x = Ax + Bu 


y = Cx 


8 



where 





(16) 


(17) 


and 

C = [u 3 [ Oj <2^! 0] 

with the NxN matrix defined by 

m-l a (Uy-ej)” 1 


(18) 


(19) 


and the 3x3 matrix defined by 

M 2 - (U3-J6)' 1 (20) 

Note that for M r <f» — U N the matrices <$J and are symmetric, although JS and M 2 
are not. The inversion in Eqs. (16) and (17) has been accomplished by means of 
the partitioning formula in Eq. (191) of [2]. 
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As proven in Appendix B, one can interpret J and" & in terms of primitive 
definitions and obtain the physical interpretation 



where 1° is the inertia matrix of the primary body referred to its own mass cen- 
ter. This interpretation guarantees that is nonsingular, which as shown in 
Appendix B implies that is nonsingular. 

CONTROLLABILITY 


The system characterized by Eq. (14) is said to be completely controllable 

[6] if any initial state x(t Q ) can be brought to any finite state x(T) in a 

finite time interval T-t by some control function u(t). 

A basic controllability theorem [6] indicates that the system of Eq. (14) 

is completely controllable if and only if the controllability matrix, 0 , as 

c 

defined by 


Q c = 


pi . -rt i ■ 2 i i A 2N+5 
B i AB i A B i ... j A B 


( 22 ) 


ii it 

has full rank, which means rank equal to the dimension of x, namely 2N+6. (Note 
that % haS dimensions (2N+6) x (2N+6)3). With the notation 


B 2 = m 2 +JM 1^ 
B 4 A 6M 2 +M 1 0 


C 

T * 


I 

% 

i 

I 


we have 


’ 0 i B 2 i 

o ! 

i 

J(-M 1 a 2 )B 4 j ... 

i 2 N+2 

jJ(-V ) b 4 

B. ! 0 ! 

2 i i 

J(-M 1 a 2 )B 4 | 

0 j ... 

1 

1 0 

i 

0 ' 
i 

i — 

' , M rt 2.N+2 
j (-MjO ) 4 

_ B 4 i 0 i 

(-Mi0 2 )B 4 I 

0 j ••• 

j 

! o 


(23) 
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The rank of Q c is twice the rank of the (N+3)x3 (N+3) matrix Q', as given by 


Q' = 
c 


B 2 j J (-M 1 a 2 )B 4 I j (-Kj^a 2 ) 2 B 4 j . . . j J (-M 1 a 2 ) N+2 B 4 


_B 4 i C-M 1 a 2 )B 4 j (-M^a 2 ) 2 b 4 I ...j C- Ml a 2 ) N+2 B 4 


(24) 


In evaluating the rank of it is convenient to premultiply by a matrix V 
which has j V } =1 and which results in a block triangular form for Q^V. By 
proceeding in three stages one can determine that the desired matrix is 


V = 


U 3 j 0 


u 3 i -JM‘ 


° 3 ! 0 1 

.T -1 *! 

-+ cV 1 i d n 


Tl 

i 

i 1 


1 

1 d*. 

1 

! L_ i 


(25) 


Since the determinant of a product is the product of determinants, the require- 
ment | V | = 1 is satisfied, and premultiplication by V won’t change the rank of 
Q\ Moreover, the product VQ^ becomes (with the definitions of B 2 , B 4 , and 


VQ’ = 
c 


I o 

i 


-<J) T m ^ 1*!' u 
* c 2 1 J U N 


u 3 ; -jMj 


0 j u 

j N 


B 2 j j(-M 1 a 2 )B 4 j ... j J (-m^ 2 N+2 B 1 


^ c 1 * j “l 1 ~ M i a2)B A | • ■ • I M^ 1 (-M 1 a 2 ) N+2 E 


u. 


I o 


T 1 A 

-+ W 1 1 °NJ 


1 

k 

O 

• 

* 

•# 

o 

1*1 

L c i 

M~ 1( -Mia 2 )B 4 | ... j M“ 1 (-M 1 a 2 ) N+2 B 4 


M. 


T * -1 * 

<t> <? c (I -m 2 a i m 2 ) 


0 




where 

^ [ B 4 1 (-V 2)B 4j ( -«l <l2 >S! 




(26) 


(27) 


Eq. (21) guarantees that M 2 is nonsingular, so that the first three columns of 
VQ' c In Eq. (26) are independent, and by virtue of the null partition these three 
columns are independent of all others in VQ\ Moreover, the possibility of 


11 



X 2 

normalizing <j> so as to make a the diagonal matrix o% squares of nonzero 

2 

natural frequencies of the appendage indicates that 0 is nonsingular. Thus 
the controllability condition becomes 

Rank(Q’ N ) = N (28) 

In the interpretation of this result, it becomes convenient to rewrite 
from Eq. (27) in terms of 

A 6+<J> T j5? I* (29) 

The identities [4] 

M 2 = (Ua-J 5 )" 1 = U 3 +J(U-6J)~ 1 (5 = l^+JM.,6 (30) 

and 

M x = (U^dJ)' 1 = = U n +6M 2 J (31) 

permit recognition of the further identities 

JM X = J(U n +6M 2 J) = (U 3 +J6M 2 )J = (M” 1 +J6)M 2 J = M 2 J (32) 

and 

« 6+6 JM^ 6 = 6-¥6K^}6 - M^6 (33) 

and 

M x = U n +(5M 2 J = U n +M x 6J (34a) 

and 

\ = U n +6M 2 J = U n + 6JM 1 (34b) 

and these permit to be written as 

B 4 = = M x 6 c (35) 

and thus permit to be written in terms of 6 t A further simplification is 
afforded by multiplying by the block diagonal matrix with blocks 
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> (-U^)^, ... which removes the explicit minus signs from Eq. 

(27) without changing its rank. If this product is called Q , then the neces- 
sary and sufficient condition for complete controllability becomes 


Rank(Q CN ) = N 


where 


*CN 


- I M. 6 ! 

1 Cl 


<V 


j( Ml a 2 ) H+ \« ( 


(36) 


(37) 


Before attempting to interpret this criterion in useful physical terms, we 
establish a parallel criterion for observability. 

OBSERVABILITY 

A state x^(t) characterized by Eqs. (14) and (15) is said to be observable 
[10] if knowledge of u(t) and y(t) over a finite time segment t Q <T£t completely 
determines x^ (t) ; if all states in a system are observable the system is said 
to be completely observable. 

A basic observability theorem [8] indicates that the system of Eqs. (14) 
and (15) is completely observable if and only if the observability matrix Q q , 
as defined by 


q o .[c t ! aVI <a t ) 2 c t ] j (A T ) 2N+5 C T 


(38) 


has full rank, which means rank equal to the dimension of x, namely 2N+6. 
Substitution from Eqs. (16) and (18) , with the notation 

6 A <(> T M , ii>5+<]>'V I* = I* 

0 = 0 0 


(39) 


(to be compared with Eq. (30), noting the eigenvector normalization - U^) , 


so that, with the definition of J, 
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J T +(^ T <?>) T = (t> T M r cf>) T <S<I*“ 1 ) T -Nf) ,i \sf 7 - ((f) r M , 4)6-hf> T ^ , ( 


*-l 


*-l 


T 

produces for the observability matrix (in which = M^) 


1 

0 

1 1 
: o ; 

0 

\ 

1 

0 

1 

J 

) 

0 

0 

U 3 

> i 

! 0 j 

0 

1 

r 

0 

r 

I 

l 

0 

T » 

0 r o i 

0 

i 2 » 

| <-o *0(5 I L ) j 

1 1 o 1 

0 

J 

J 

f 

(-a 2 M 1 ) N+2 (6 o I*"^) 

l 

j 

i 

0 

r 

0 ! 

i 

.T , 

* -''o 

j o ! <- 

! I 

c 2 M, ) (<S I* -1 ) 
I o 

i 

t 

1 

i 

f 

0 

|'(-o 2 m 1 ) 

N+2 (6 I*" 1 ) 
o 


The rank of Q is twice the rank of O', where 
o o 


( 40 ) 


U_ 


i* = 


0 


' (-a 2 M ) (6 I*" 1 ) ' (-a 2 M ) 2 (<S I*” 1 ) 
O 10 1 o 


\ 

' (-a 2 M )^ 2 (6 I* - *) 

1 O 


(41) 


and since the first three columns of are independent of each other and of 
all other columns the rank of is three plus the rank of the residual matrix 
obtained by deleting these three columns and the three top empty rows from 
^o* ® ecause jl J 4 0, and [c 2 j ^ 0, this residual matrix has the same rank 
as the matrix 


Q 0N A £Vo 1 ! ••• 1 (42) 

(Note that, as in the controllability case, minus signs have been removed, 
without influencing rank.) 

Thus the necessary and sufficient condition for complete observability 
becomes 

Rank(C W = N (43) 

The observability conditions in Eqs. (42) and (43) are now very similar 

to the controllability conditions in Eqs. (36) and (37), differing only in the 

exchange of 6 for 6 . 

o c 
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INTERPRETATION OF RESULTS 


The reduction of the controllability matrix to (Eq. (37)) and the obser 
vabllity matrix to Q qn (Eq. (42)) implies that controllability and observability 
of the system of dimension 2N+6 in Eqs. (14) and (15) can be inferred from that 
of the system of dimension N given by 

v * MjCT^v+iij^u (44) 

\ 

w = d^M-v (45) 

o 1 

With this interpretation comes the realization that simpler necessary and suf- 
ficient conditions for controllability and observability are given respectively 
by the requirements for full rank (N) of the matrices 

Q C N ^ tVc ! (M l a2) Vc j J <\ 0 2 )^\ 6 C ) (46) 

and 

i ‘Vo! V’V.,!- ! (47> 

These conditions and the relationships between them are most interpretable in 
physical terms when special cases are considered. 

When all of the sensors and actuators are attached to a single rigid body 
in the model, we can call that body the primary body, so that the location 
matrices 3 ? and & are both zero, 6=6 = 6, and the controllability and 

observability conditions become identical, and simplify greatly. In this case 
Q cn ln E( i* C 46 ) and Q qn in Eq. (47) can be written with Eq. (33) in the form 

^CN = ^ON * i ^ M i a *** i ^ ^2^ (48) 

Since by Eq. (21) the common factor M 2 is nonsingular, the rank of the matrix 
in Eq. (48) is unchanged by removing M 2 throughout ; thus the controllability and 
observability criterion becomes 
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(49) 


Rank [ 6 j (MjAfS j ... | (M ^ 2 )* 1 1 $] = N 

Eq. (A9) simplifies dramatically when subjected to the following corollary 
of Kalman's canonical structure theorem [7]: "If a system characterized by the 

pair (A,B) is completely controllable and F is any matrix of appropriate dimen- 
sion, then the system characterized by (A+BF, B) is completely controllable." 
(This corollary is not proven in reference [7], so a proof is included here as 
Appendix C . ) 

Eq. (34b) permits Eq. (49) to be written as 


Rank [ 6 j (cAsjM^ct 2 ^ | ... j (c^+fiJM^ 2 ) 15 1 6] = N 
or, with the noted corollary. 

Rank [ 6 J o 2 6 j ... } (a 2 ) N * 1 fij = N 

A typical Nx3 partition in the preceding matrix is 


(50) 


2U 

a o = 


2i 


.21 


a. 


2i 

N 


JJ 


2i x i ! 

21.1 

! 2i .1 

a i 6 i ! 

a l 5 2 

i a i 6 3 

_2i .2 I 

°2 6 l! 

« 1 

-21..2 
°2 6 2 

■ 

|°2 6 3 

1 • 

• I 

. 1 

rr 21 ^ < 

ff H 5 i! 

• 

« 

°H 6 2 

1 i 

J * 

! 21, N 

: °n s 3 


The Nxl columns in Eq. (50) can be interchanged without changing the rank, so 
that Eq. (50) becomes 


Rank 


2.1 

* 1*1 

2.2 

a 2 6 l 


S 1 Vi 


„2<N-1K2 


2(N-1) .N 
a N 6 1 


.N 


°1 fi 2 

2 (N-l) .2 
CT 2 6 2 


2(N-1)N 
°N 5 2 


a — 2 (N-i) .i 

j 0 1 6 3 

.2 .... 2 ( N - 1) .2 

« 3 o 2 « 3 


,N 2(H-1).N 

6 3 °N S 3 


=N (51) 
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where 


a o o 

Rank[A 1 a j j A 3 $] = Rank[A., | A, j AJ 0 o 0 

ooa 


A t A 


1 J. 4 

1 a N a N 


2(N-1) 

1 

2 (N-l) 


2 (N-l) 


Since the determinant of the Vandermonde matrix a is given by [8] 


lal = H(a 2 -a 2 ) 


1 £ i < j £ N 


then 3 is nonsingular if for all i + Oy Under this supposition Eq. 

(50) is equivalent to 


Rank[A 1 j A 2 j A 3 J 


Rank[A j A j A j 


= Rank 


6 6 0 


0 «V 


thus for for all i ^ j the condition 


6 6 '+ 0 


i “ 1, ... , N 


is necessary and sufficient for both controllability and observability. 


17 



In the case of repeated roots, we may assume that the first p roots are 
repeated and the remainder all distinct, and then re-examine Eq. (51). Now 
the first p rows are independent if and only if 6 1 , ... , are independent. 
(Note that this is impossible for p>3, or in general for p exceeding the dimen- 
sion of 0, which establishes the number of control axes in the problem.) If 
the top p rows are deleted from Eq. (51), then what remains has a rank at least 
as large as that portion of it given by a new version of Eq. (52) with N re- 
placed by N-p and A and o beginning with 0 and S? +1 rather than a, and 6 1 . 

1 P+1 i 1 i 

Thus in the case of p repeated roots a. , ... , a it is sufficient for control- 

1 P T 

lability and observability that 6 1 , ... , 6 P be independent and 6 i 6 i j£ 0 for 
i = p+i, . .., N. (The authors suspect that this is also a necessary condition, 
but the proof has eluded them.) 

Eq. (53) admits an appealing physical interpretation when written as 


tr(6 i 6 1 ) £ 0 (54) 

T i T i 

since 5 5= £ 5 6 has for the normalization «= previously been in- 

teirpreted ([2], pp. 69-70) as the inertia matrix of the appendage referred to 

the system mass center c plus the difference in the inertia matrices of the 

>T 

primary body referred to c and to its own mass center. The quantity 6 1 <$* has 
been called the "effective inertia matrix" of the 1 th mode, and has been recog- 
nized ([2], pp. 69-70) as a measure of the coupling between i th mode vibration 

. .T 

and primary body rotation. If <$V =0, the i th mode would for output 0 be de- 

leted from the model by coordinate truncation; thus the truncated system is 
completely controllable and observable. 

This success in dealing with the special case in which sensors and actuators 
are all attached to one rigid body (here called the primary body) has not been 
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matched for the general case of arbitrarily distributed sensors and actuators. 
Since the observability condition (Eq. (38)) depends only upon C and A while 
the controllability condition (Eq. (22)) depends only upon B and A, one could 
cope with the special case in which all sensors are on one rigid body and all 
actuators on another by applying Eq. (53) twice, with two different selections 
for the primary body and two different interpretations of all of the symbols, 
but this case is less interesting than the common primary body case solved here 
or the general case in which sensors and actuators are both distributed over 
the body (as in the case of the telescopic sensors on the NASA Large Space 
Telescope) . For this general case the controllability and observability con- 
ditions are full rank requirements for the matrices in Eqs. (46) and (47) res- 
pectively. These requirements can be placed in slightly simpler form by observ- 
ing that |M^[ ^ 0 (Appendix B) and factoring to the left; necessary and suf- 
ficient conditions for controllability and observability then become respectively 


and 


Rank I 


Rank 


[ 6 c i (0 S )6 ci--- i 

£.1 <‘Vol- 

These criteria correspond to the state and observation equations 


N 


= N 


v = a 


W = 0 V 

o 


2 M 1 v+< 5 c U 


(55) 


(56) 


(57) 

(58) 


For purposes of preliminary analysis, it may be useful to consider a flex- 
ible appendage model characterized by a single mode. Eqs. (55) and (56) then 

require simply that the 1x3 matrices and S q be of rank one, or the scalar 
equivalents 
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(59) 


and 


6 6 T 

c c 


0 


n T ^o 

o o 


(60) 


Experience with specific examples (see Appendix A) and with the special 
case represented by Eq. (53) suggests that when many appendage modes are in- 
cluded the satisfaction of Eqs. (59) and (60) for each mode is necessary but 
not sufficient for complete controllability and observability respectively, 

with an additional frequency criterion stipulating further criteria in the 

2 

event that eigenvalues of O are repeated. This proposition has however 
not been proven; the corollary which produced Eq. (50) from Eq. (49) and sim- 
plified the special case governed by these equations is not helpful in the 
general case. 

CONCLUSIONS 

Necessary and sufficient conditions for controllability and observability 

have been established in the form of full rank requirements for Nx3N matrices 

for systems characterized by state equations of dimension 2N+6 . For the 

special case in which all sensors and actuators are attached to the same "pri- 

mary" rigid body, the controllability and observability conditions are identical 

. .T 

and reduce to the requirements that N scalars 6 1 6' L (i=l, ... , N) be nonzero and 
either appendage frequencies cr^, ... , cr^ be distinct or effective inertia 
matrices 6 1 6 1 be independent. 

With these conclusions come a new reason for the study of controllability 
and observability. This investigation was initially stimulated by the realiza- 
tion that much of the progress in optimal control theory depends upon the assump- 
tions of complete controllability and observability [9]; if we are to address 
the problem of flexible spacecraft control with the techniques of modern control 
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theory we must first determine whether or not our system is completely control- 
lable and observable. It has been apparent from the outset that this issue is 
clouded by the realization that in the flexible spacecraft modeling process 
there is an ambiguity stemming from the foreknowledge that not all possible 
modes of vibration are dynamically significant or even computationally available, 
and many must be truncated from the system description prior to simulation. The 
selection of which modes to truncate and which to preserve is a major decision 
in flexible spacecraft simulation. Now we can see that the preceding conditions 
for controllability and observability not only provide answers required for the 
application of the techniques of optimal control, they also provide a new ration- 
ale for modal coordinate truncation, as required for any system simulation. 

It has been recognized previously ([2], pp. 69^70) that the "effective 
inertia matrix" of each mode affords a rationale for truncation, since this 
matrix provides a quantitative measure of the dynamic coupling between vibra- 
tion in that mode and primary body rotation. Now we see from Eq. (53) that 
the trace of this matrix is a formal measure of controllability and observabil- 
ity when all sensors and actuators are on the primary body. Moreover, we find 
that in the more general case in which controllers and observers are distributed 
over the body, we have a new truncation criterion in the controllability and 
observability conditions available from Eqs. (55) and (56). As we entertain 
various truncations, our objective is to obtain a system of minimal dimension 
consistent with the preservation of the fidelity of the output matrix 

Z “ DX (61) 
which describes the salient aspects of the spacecraft mission performance. 

Since the control law for u in Eq. (11) depends upon the sensor observations 
y in Eq. (13), it is necessary that all significantly observable states be re- 
tained in the coordinate truncation process. Although in general it is possible 
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for an uncontrollable state to influence z and hence be among the states which 
must be retained in truncation, in the special case of states describing flexible 
body inodes experience suggests that only controllable states will influence z. 
Under this assumption, it becomes desirable to truncate all those uncontrollable 
states which are also unobservable. This leads to the adoption of the design 
objective of configuring sensors so that any uncontrollable states which do not 
influence z are also unobservable; this strategy of uncontrollable mode isolation 
permits in many flexible spacecraft control applications the truncation of all 
uncontrollable states, which are by design also unobservable. The result is 
often a completely controllable and observable system model, which is more amena- 
ble to formal optimal control, simulation, and conventional control system design. 
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APPENDIX A. EXAMPLES OF CONTROLLABILITY AND OBSERVABILITY CALCULATIONS 

Example 1. Consider single axis rotation of a spacecraft with symmetric append 
ages of continuous thin beams as shown in Fig. Al. 



The rotational motion of the model is characterized by Eqs. (11) and (12) with 
the observation equation (13). However, the assumptions that the rotational 
motion is limited to the single axis and that the mass distribution is contin- 
uous (so that each mass element has infinitesimally small moment of inertia) 

necessitate a slight modification of definitions of 6, J, <(), SB and SB • Fur- 

c o 

thermore, mode shapes are assumed to be given by a function of £, i.e., <j>%) 
which is at least once differentiable with respect , to 5 over the intervals 
d<[C|<L+d, or more specifically, 

0 ,0<|C|<d 

< »d<£<L+d (A-l) 

t - (L+d) <£<-d 

Preceding page blank 
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<D j (0 = 



Under these assumptions, we may express (S' 3 (scalar) in terms of mode shape and 
inertial parameters as follows. 


fiJ - - 


m 

L 



<J>%)£d£ 





(A-2) 


For the purpose of illustration, we further assume that only two modes are 
retained (N=2) . The first modal coordinate, is associated with an asymmetric 
mode which satisfies 


th-O - 

and the second is with a symmetric mode which satisfies 

4 l(-£) = 4> 2 (£) 

The modal deformation coordinates thus defined are illustrated in Fig. A. 2. 




Figure A 2. Asymmetric and Symmetric Modes. 


Corresponding to these modes, 6 J ’s (j=l,2) are given by 


/ L+d 


(0£d£ 


6 2 = 0 


We examine the following three cases. 


(A- 3a) 
(A-3b) 
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a) If both sensor and actuator are located on the primary body, the rotational 
equation is characterized by 
.. *1 .. 


f e - ^ n = u 

i 

• • ^ 1 mm 

n i i n i = 0 


V 0 ^ 


(A-4) 

(A-5) 

(A-6) 


with observation equation 

y - 9 (a-7) 

Since Eq. (A-6) is independent of u, 0 and and does not affect Eq. 

(A-7), we can readily identify rj^ as uncontrollable and unobservable. 

Thus, we may limit our analysis to the truncated system of equations in which 
Eq. (A-6) is deleted. 

Direct approach: For comparison, we examine the controllability and 

observability matrices directly. In terms of the state variable 
x A [9,e,n 1 ,ri 1 l the state equation is x * Ax+Bu with 


A = 


0 

1 

0 

0 


0 



*1„2 
-6 o 




0 

n 

1 

A 


1 

y 

iV) 2 

U 

and B = 

I-C6 1 ) 2 

0 

0 

0 

l 


0 

0 

n 

* 2 
- I0 i 

A 


6 1 

U 

i-W 1 ) 2 

u 


Li-(6VJ 


and with the observation equation y = Cx with 
C = [1 0 0 0] 

Now, the controllability matrix is 
Q c = [B | AB | A 2 B j A 3 B] 
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and the observability matrix is 


Q = [C T | a t c t ! (a t ) 2 c t | (a t ) 3 c t ]. 


* 12 

Matrix algebra produces, with a A (I-(6 ) ) (>0) , 


a Q = 


and 


aQ = 
x o 


0 

a 


.1 

6 a 


a 

0 

0 

0 


a 


0 6 a 

0 


0 

a 

0 

0 


0 


T *K 2 
-I 6 ff 1 


0 

0 

- 6^1 


,. 1*2 2 

-(o ) 

o 

*.i 2 

-I o 


0 

0 

0 


o 


Consequently, 

|Q | » (-6 1 ) 3 af/a 2 

C 1 

|Q | = (fih^/a 2 

' x o 1 1 


Since a / 0, a>0, the full rank conditions reduce to 
S 1 # 0 

for controllability and observability. 

The result (Eq. (53)) on the controllability and the observability 

produces the same conclusion as follows. 

2 

The fact that 6 =0 implies that the state including T) ^ is not com- 

pletely controllable, and is not completely observable either because the 
sensor is attached to the same body as the actuator. 

If n 2 is deleted from the state, then Eq. (53) simply requires 
6 1 ^ 0, noting that N = 1, 
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b) If the actuator is attached to the primary body and the sensor is attached 
at £ * K as shown in Fig. Al, then it reads 


1 ^ Of 

y = 8 +[<fr <r ] 


(A-8) 


where for simplicity of notation we write 




S * 5. 


k <*>!'» j = 1, 2 ' 


In this case, Eqs. (A-4)-(A-6) remain unchanged so that the, control- 
lability condition is the same as the previous result. 

The observability is to be discussed with respect to the state in- 
cluding n 2 in view of the difference between Eqs. (A-7) and (A-8). By 
definitions of Eqs. (19) and (39), 


M, 


- 




(?) 


-1 


I 

a 


0 

1 


L 

1 


L J 

J 

L 

= 

H 

J 

+ 

.2' 

H 

* 

l * l r 
6 X +I <j) x 
s 

* 2 ’ 


L° J 






and consequently, 


M .6 = 

1 o 


* 


* 

J s 1 

a o 

» M 1 a 2 = 

K o 

iV 

o 


_° «?. 
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Thus, we can construct defined by Eq. (47) 


<0N 


■ * 

a o 

* .2 

I 

o 


* 2 

° 2 A 


Since is a square matrix, the full rank property is simply examined 
by its determinant, i.e., 


*2 * 

|q | = - — (S 1 S 2 (a^ - ^-oh 

1X1 a oo x 2 a 1 ' 


Thus, observability requires three conditions: 

(i) S 1 = 6 1 +I% 1 V 0 

o s 

? * 9 f 

(ii) Q = I (p ^ 0 

o T s 

* 

(iii) a 2 - ft 0 


We see that the conditions (i) and (ii) are exactly what Eq. (60) requires. 

2 * 

Especially, the second requirement which reduces to 4> ^ 0 is physically 

s 

2 1 

clear because if <J> a =0 then Eq. (A-8) is not affected by n 2 as in the 
previous case. 

The condition (iii) deserves particular attention. This prohibits the 

2 9 

diagonal elements of M^a to be equal or the eigenvalues of M^a to be 

repeated* Eq* (44) should be referred to with consideration that in case 
2 

of diagonal t repeated eigenvalues are not accepted for the single input 
to be either completely controllable or observable# (For similar 
discussion, see Ref . [ 9] ) * 
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c) Assume that the sensor is attached to the primary body and a pair of gas 
jets are attached on the beams at £ = l and £ = - Further, assume that 
two actuators synchronize but their forces are not necessarily of equal 
magnitude, i.e., the one at £ = - £ produces a force f^ which is given by 
f£ = c'f^ with - 1 < c '<1 where f^ is a force produced by the jet at E, - &. 
Such dependency comes from the minimal actuator model. Then Eqs. (A-4) - 
(A-6) are rewritten as 



■ —^(l-c')f = u 

I 


r)-HJ n-<50 


where c 



— — 


~ — 

<{£(*) (1-C) 

f, - 



1 


^(JD(l+c*) 


^(Jl)l*c/£ 


By the definition of Eq. (29), 



6 1 

c 


6 1 +(J>^(J2) l*SL 

6 = 


= 


c 

o 


2 * 


<$ 

c 


<j>+ a n 


Similarly as previously, we can construct Q 

"A 1 (i)w 


defined by Eq. (47) 


CN 


* „2 
I 6 

c 


*2 2 

I o,o 

1 c 


and 





a 


5 c 6 c (6 2- 
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Thus complete controllability requires three conditions: 

(i) s* = t o 

(ii) & 2 c = <**(*) I V* * 0 

(iii) a\ - f o 

The condition (iii) is identical to that of (b) . (i) and (ii) are what 

Eq. (59) requires. The condition (i) together with (A3-a) implies that 
if the mode iJ^CO satisfies the integral equation 

/•L+d * 

- ^(OSdS+j <£(i) - 0 

then is uncontrollable. (Equivalently, if the actuator location 
parameter, satisfies this, then t\ is uncontrollable.) 

The condition (ii) requires 

4- 0 and c^O 

c = 0 implies c T = -1 and ^ which is realized by a pair of equal 

magnitude jets. From the practical point of view, it is rather desirable - 

2 

to locate the jets at Z satisfying “0, because it can eliminate 

the effect of which causes an unnecessary vibration due to some im- 
perfection of the jets (c ^ -1 or f^ t -f^) . If this is possible, then 
we can truncate for this particular location, although otherwise we 
must retain it. 
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APPENDIX B. INTERPRETATION OF 
By the definitions following Eq. (13) , 

U 3 -J6 = U 3 -I* S T (4> T M'<f>)S * I* [ ($ T M t tf>) 6 ] 

Moreover, in view of the symmetry of M f and the definition of <5, 
<$V M’4>6 = A T <t(/M , (j))" :L (4> T M , cf l )(0 T M , ^) _1 4) T A 

= A T ^ -1 M ,_1 4) T <J) T A ^ a t m ,_1 a 


(B-l) 


(B-2) 


From [2], Eq. (277), with minor notational revision, 

M ' A "(V Voo^ > ®-3> 

and 

A 4 ^oiT^UO 5 (B-4) 

where M is the block diagonal 6nx6n matrix containing in alternating sequence 
along the main diagonal the mass and inertia matrices of the nodal bodies, as 

designated by the 3x3 matrices m , I , m , I , , m 11 , I n . The symbols 

denotes the total vehicle mass, and the 6nx3 matrix binary operators and 
2 0U are defined in terms of 3x3 zeros and unit matrices by 

z uo ■ °! U 3 j OJ ... ( U 3 ! 0] T S (B— 5) 

and 

J.OU = *-° I u 3 i 0 ! U 3i • * • i 0 ! U 3^ T (B-6) 

Finally, the ft in Eq. (B-4) is defined in terms of 3x3 zeros and the matrices 
in Eq. (5) by 
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R A 


R 

0 

0 


0 0 


0 0 
0 R 2 


R n 0 


0 0 


(B-7) 


Thus Eq. (B-2) becomes, in view of the symmetry of M, 




i T M* 1 A = 
with the identity [4] 


‘"sn-Vuo'^'”' 1 ■ U 6n + 1 V '^lo^ 


and the recognition that 


n 


n 


E UO ME UO ■ , E / - , E "A 
1=1 1=1 

where nu is the mass of nodal body i, we can identify the scalar 
n -1 

(1 - S m Mt) A = e ^T/m 
i-1 

where m Q is the mass of the primary body, and write (noting the skew-symmetry 
of R) 


- (Z* u+ zJ 0 R)M(u 6n+ ^VCVKdo, 


(B— 8) 


With the identities 


z o^ l m - 0 


E So ME ou ’ 0 
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^uo^ou “ 0 


- 0 


^OU '. E 1 lJ 
J=1 


3=1 J 


Z UO fi ^ iE UO " E UO MftS UO 


Eq. (B-8) becomes 


n 


A T M ,_1 A -21 j - E m 4 J[ j R j - [E m J^] 2 /m 


j=l 3=1 j 


3-1 j 


(B-9) 


The first two summations can be identified by the inertia matrix reference point 
transfer theorem [5] as the inertia matrix of the set of n appendage nodal bodies 
referred to the system mass center c. Thus, if I' is the inertia matrix of the 
primary body referred to c, then, with Eqs. (B-9) and B-2), Eq. (B-l) becomes 

.-1 n 


U,-J<$ - i" [I* + (E m.S j ) 2 /m ] 
3 j=l J 


(B— 10) 


The sum Em^R^ is by system mass center definition equal to times the vector 
from c to the mass center of the primary body; thus by the noted reference 
point transfer theorem the expression in square brackets in Eq. (10) is the 
inertia matrix 1° for the primary body referred to its own mass center. This 
matrix is symmetric and positive definite, which guarantees that the matrix 
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(B-ll) 


* o 

U 3 -J6 = I I 

is nonsingular, and the matrix 

-1 o -1 * 

M 2 = (U 3 -J6) 1 - 1° I 


(B-12) 


exists, so that 


| U 3 -J6 i + 0 


(B-13) 


From the relationship [4] 


! ? i 

F 11 ! 

| = det — 1 
_ P 21 1 

_ P 12_ 

P 22 


l p nl 

1 1 P 22 _P 21 P 11 P 

12 ^ = i P 22 1 

l P ll" P 12 P 22 P 21 

with P n = U 3 , 

P 22 = U N* P 12 

= J and = 

<$, we see that 


|U N - J <$| = |u 3 -6j| 


(B-14) 


(B-15) 


so that the nonsingularity of is assured by the nonsingularity of 
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APPENDIX C. PROOF OF CONTROLLABILITY COROLLARY 


Corollary ; If a system characterized by the pair (A,B) is completely state 

controllable and F is any matrix of appropriate dimension, then the system 

characterized by (A+BF,B) is completely state controllable. 

Note : This corollary is applicable whenever there exists a constant linear 

feedback law, such that the control variable u(t) can be written in the form 

u(t) = u (t)+Fx(t) 
c 

where F Is a constant matrix. Then the corollary states that the system con- 
trollability is unaffected by the linear feedback term Fx(t). 

Proof 

Define a tridiagonal matrix V by 



Q c A [B i ( A+BF) B j (A+BF) 2 B j ! (A+BF) n_1 B] 

so that 

VQ c = [B j AB j A 2 Bj — | A n_1 B] 

Since |v| =1, the rank of is the rank of VQ^, and the corollary is proven. 
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CHAPTER II 


Commutativity of Coordinate Truncation and Transformation 

•f* 

Matrix Inversion for Flexible Spacecraft Dynamic Analysis 


P. Likins* 
and 

Y . Ohkami 


ABSTRACT . Necessary and sufficient conditions are established for the com- 
mutativity of matrix Inversion involved in certain coordinate transformations 
and the coordinate truncations that are essential to practical structural 
dynamics. Computational constraints demand that truncations precede the 
generation of explicit transformed equations, while mathematical arguments 
imply the opposite sequence. Results jeopardize the utility of two of the 
three transformation procedures considered here for rotating flexible bodies. 
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introduction 


It is now commonplace to characterize spacecraft with flexible appenda- 
ges in terms of hybrid coordinates, consisting of discrete coordinates for 
rigid bodies in the system and distributed or modal coordinates for the 
flexible appendages* The selection of mode shapes and frequencies associ— 
ated with the modal coordinates is usually accomplished by adopting a finite 
element model of the appendage and writing equations of motion for the small 
vibrations of the appendage with respect to a base which has a nominally 
constant (perhaps zero) inertial angular velocity; the eigenvectors of 
these equations then provide the mode shapes and the eigenvalues provide 
the natural frequencies of appendage vibration. This information is used 
in various ways to construct a transformation matrix which transforms the 
discrete nodal coordinates of the finite element model into distributed 
coordinates; these then are reduced in number by a process of coordinate 
truncation, which sacrifices mathematical rigor for computational feasibility 
(hopefully without doing violence to the salient features of the mathematical 
model). This paper addresses an aspect of the truncation problem which 
arises in the case of modal analysis of an elastic appendage on a rotating 
base. Several alternative transformations have been proposed for this case, 
but comparison among them is difficult without detailed numerical studies. 

In this paper we consider three alternative transformations, two of which 
involve the inversion of the transformation matrix. Whereas mathematical 
arguments demand that this matrix be inverted prior to truncation, practical 
considerations demand that truncation precede inversion (which then requires 
a pseudo— inverse) , In this paper we establish necessary and sufficient 
conditions for the commutativity of truncation and inversion, in order to 

Preceding page blank 
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define the formal limits of these two procedures. We conclude with recom- 
mendations for the selection of coordinate transformations for elastic 
appendages on rotating bodies. 
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EQUATIONS OF APPENDAGE VIBRATION 

As shown in Eq. (64) of [1], the equations of vibration of a finite 
element model of a flexible appendage on a rotating base have the structure 

M'q + D*q + G*q + K*q + A*q = L* (1) 

where for an appendage with n nodes q is the 6nxl matrix of small deforma- 

i » t 

tions relative to a nominal (deformed) state, the matrices M , D , and K 

I f I 

are symmetric, and the matrices G and A are skew-symmetric, with M positive 
definite and hence nonsingular. These equations are applicable whether the 
appendage mass is concentrated into rigid nodal bodies (as in [2], Eq. (140)), 
or also distributed over the finite elements of the model (as in [1]). The 

i 

matrix L in Eq. (1) holds the vibration forcing functions, which include 
base motions, but the coefficient matrices in Eq. (1) also depend on the 
inertial angular velocity of the base. If and only if this velocity 

experiences only small deviations from a nominal constant can these coef- 

\ 

ficient matrices be treated (after linearization) as constants. 

COORDINATE TRANSFORMATIONS AND TRUNCATION 
FOR NONROTATING ELASTIC APPENDAGES 

It is well known [3] that when the base is nonrotating and the appendage 
is undamped the homogeneous counterpart to Eq. (1) becomes 

i i 

M q + K q = 0 (2) 

and that the transformation to modal coordinates 

q - <j>U (3) 

T 

followed by premultiplication by <f> produces an uncoupled system of scalar 
equations 

2 

n ± + CK Uf = ° i = 1, ...,6n (4) 

2 

where the columns of $ are the eigenvectors of Eq. (2) and represents the 
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til 

i eigenvalue. The inhomogeneous counterpart to Eq. (4) may be written as the 
scalar equation 

n i + «J 1 2 n 1 = (<f> T L , ) i i = 1 6n (5) 

T * th 

where (4> L ) is the i element of the 6nxl matrix in parentheses. The 
vibration equations in Eq. (5) could theoretically be used in conjunction with 
the vehicle rotation equations to fully characterize the spacecraft attitude 
dynamics, but in practice this is not done because most of the modal coordi- 
nates (i = l,...,6n) have no significant influence on the motion of the 
primary rigid body to which the flexible appendage is attached. Thus the 
modal vibration equations are ignored except for a small number N < 6n 
(typically N ranges from one to twenty and 6n is measured in hundreds or 
thousands) and Eqs. (3) and (5) are replaced by 

q - (6) 

and 

n ± + C^ 2 !^ = (4> T L ) i i = 1,...,N (7a) 

or 

r\ + an = 4> L (7b) 

where n is the Nxl matrix of truncated modal coordinates, a is the NXN 
diagonal matrix of corresponding squared natural frequencies, and | is a 
Nx6n matrix whose columns are the eigenvectors corresponding to those modal 
coordinates which are retained. In practice one computes only the N eigen- 
vectors in 4> and the corresponding eigenvalues, which normally include the 
lower eigenvalues; calculation of all 6n eigenvalues and eigenvectors would 
be not only prohibitively expensive but also computationally infeasible for 
typical values of n. 

Thus in the special case of Eq. (1) represented by Eq. (2) the trans- 
formation matrix <p (in Eq. (3)) need not be inverted; since the required 
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transposition of <}> is commutative with truncation, there is no mathematical 
obstacle to the pragmatically motivated practice of truncating <{> (as in Eq. 
(6)) before transposing it (as in Eq. (7)). 

GENERAL PROBLEM OF TRANSFORMATION 
AND TRUNCATION OF COORDINATES 

Although other special cases of Eq. (1) permit the use of special trans- 
formation matrices which circumvent all matrix inversion (see [1], pp. 726-729 
and [2], pp. 46-56), we are here concerned with the more general case, for 
which there exists no transformation of the structure of Eq. (3) which 
transforms the homogeneous counterpart to Eq. (1) into uncoupled scalar 
second order equations such as those in Eq. (4). In the general case one 
must rewrite Eq. (1) as a first order matrix (state) equation before attempt- 
ing a transformation to uncoupled scalar equations. 

In terms of the state variables in the 12nxl matrix 



Eq. (1) can be written as 
+3?Q =£ 2 ? 

or alternatively as 
Q = BQ + L 

where 


I 

* ! i 


r ! » t i 

K +A 

i 0 

; m a 

O 

1 I 

> 

0 


i t i i » 

K +A [ D +G J 


( 8 ) 

(9) 

( 10 ) 


B A 


r • 

0 _____ 

-(M’)-i(K’+A’) 



U_ 

-1 
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L A 


0 

(MVV 



1 X2n 

The eigenvalues A^ 2 n and the eigenvectors n associated 

with the differential operators on Q in Eqs. (9) and (10) are the same, and 
since all quantities in these equations are real we are assured that any 
complex eigenvalues or eigenvectors appear in conjugate pairs. In what follows 
we assume that the eigenvectors are linearly independent, so that the matrix 
$ whose columns are these eigenvectors is nonsingular. The validity of this 
assumption is assured if the eigenvalues are distinct [4], and assured even 

i i 

in the presence of repeated eigenvalues if A =D =0 and the null solution of 
the homogeneous counterpart to Eq. (9) is Liapunov stable ([2], page 43). 

The coordinate transformation 


Q = $Y 


( 11 ) 


can then be used to good purpose in either Eq. (9) or Eq. (10). The first 
of these alternatives is developed in [1], with results only summarized 
here. If $ is the matrix whose columns are the eigenvectors of the equation 
adjoint to the homogeneous part of Eq. (9), namely, 

q’ + 38 T q = o 

f 

then by substituting Eq. (11) into Eq. (9) and premultiplying by $ one can 
obtain 

Y = AY + (** T j/«)"V T .2? (12) 

'T 

where ($ is diagonal and 


A = ($ , ^/$)“ 1 ($ ,T .35$) 

is the diagonal matrix of eigenvalues A^,...,A^ 2 n « 
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As noted in the context of Eqs. (5) - (7), practical considerations 
mandate the truncation of the transformed coordinates, so that the 12nxl 
matrix Y must be replaced by the 2Nxl matrix Y. It is irrelevant whether the 
truncation is Imposed on Eq. (12) or on 0 in Eq. (11) and its counterpart 0 , 
as may be established formally by partitioning 

$ 4 | ¥ 


(13) 


and 


. _» | =» 

$ 4 | $ j $ 


into the 12nx2N partitions with one overbar to be preserved in truncation and 
the remaining 12nx(12n-2N) matrices of eigenvectors to be removed in trun- 
cation, and writing 

I 

1 -1 


f I m ' 

$ 1 


¥ | I 

I 

L i J 

-i 

r -^-'T — 1 i 

¥ V 4> } 0 

|[r T . 

- M 


o ; 


($ ,T ^*)' 

\ 

- 1 l 

0 




0 

i (¥'V7 )-i’ 



(14) 

The indicated commutativity of truncation of the transformation matrix and 

*T 

inversion of $ c*/ $ is here established as a consequence of the diagonal 
structure of the latter. Obviously it is advantageous to truncate before 
inversion, because then there is no need even to calculate 7 and 1> , which 

t 

are generally of much larger dimension than $ and $ . Thus Eq. (12) can for 
practical purposes be replaced by its truncated counterpart 

Y - Ay + (?'’L/?)~W (15) 

If the transformation in Eq. (11) is introduced into Eq. (10) rather 
than into Eq, (9), then the calculation of the adjoint eigenvectors is not 
necessary, since premultiplication of the result by produces 

(16) 


Y = AY + $ -1 L 
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If truncation is imposed on Eq. (16), the result is 

Y = AY + ($ _1 )L (17) 

where ($ consists of the first 2N rows of the matrix As a practical 

matter, however, $ ^ is not generally available, since for typical values of 
n it is not feasible even to calculate all of the eigenvectors comprising 
the columns of $. In order to circumvent the problem, it was proposed in 

[2] that Eq. (17) be replaced by 

Y = W + (?) + L (18) 

where 


A ($ T $)‘' 1 0 T 


(19) 


is a pseudoinverse of the 12nx2N matrix $. This substitution amounts to the 
assumption that inversion and truncation of $ are commutative operations. 


COUNTER-EXAMPLE 


Dr. William Hooker of the Lockheed Palo Alto Research Laboratories has 
noted (in personal correspondence) that the above assumption is not always 
valid, as demonstrated by a counter-example such as the following: 

If 


then 


and 


$ = 


1 

L2 


O' 


-1 


1 

2 


3' 

4. 


-4 3 

2 -1 


whereas 


((T 1 ) = | [ -4 3 ] 


$ + = (S^)" 1 ^ = i [ l 2 ] t ($ -1 ) 
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Dr. Hooker pointed out that Eq. (18) is the formal consequence of the substi- 
tution of 

Q s $Y (20) 

instead of Eq. (11) into Eq. (10), noting that while this substitution may 
appear to be plausible its consequences are not equivalent to Eq. (17). 

PURPOSE 

Since in practical terms Eq. (17) is not directly available and Eq. (18) 
is available upon inversion of the 2Nx2N (generally complex) matrix it 

is important that we discover those conditions under which these equations 
are equivalent. This is the primary purpose of this paper. 

CONDITIONS FOR COMMUTATIVITY OF INVERSION AND TRUNCATION 
Proposition : For any nonsingular square matrix M partitioned as 



and its truncated counterpart 



the condition 


T T 
c d + a b « 0 

is necessary and sufficient for the equivalence 
M + . (M -1 ) 

where 


( 21 ) 


( 22 ) 


= (M T M)“ 1 M T 


and (M ) is the uppermost row-partition of M -1 having the dimensions of M T . 




Proof: Since [5] 


(a-bd 1 c) ^ 

-1, ,, -1. v-l 1 

-a b(d-ca b) 

-d ^c(a-bd X c) ^ 

(d-ca^b)" 1 


then 

(M -1 ) = [ (a-bd^c)" 1 | -a -1 b (d-ca~^b ) ~^ ] 

while 

H+ = (Kh I |[-M)' 1 1- 1 !- 1 ) - |-w.r Kh 

f, T . I v-l T | / T , T ..-1 T 1 

= | (a a+c c) a j (a a+c c) c J 

As necessary and sufficient conditions for (M ) = M we have 

, T .T .-IT , ~1 .-1 

(a a+c c) a = (a-bd c) 

and 

,T .1.-11 ~1 V/J -L v-l 

(a a+c c) c = -a b(d-ca b) 

The identities 


(a T a+c T c) 1 a T = (a T a+c T c)~ 1 t (a 1 )" 1 ]" 1 = [ (a T ) _1 (a T a+c T c) ] -1 


and similarly 


, T . T v-l T r , Tv-1, T T v,-l 
(a a+c c) c *[(c) (a a+c c)] 


and 


( 24 ) 


(25) 


(26) 


-a ^(d-ca ^b) * = -[ (d-ca ‘*'b) (a ^b) * 

permit the inverses of Eqs. (25) and (26) to be written respectively as 

, Tv-1, T T » , ,-l 

(a ) (a a+c c) = a-bd c 

and 

(c T ) _1 (a T a+c T c) = -(d-ca _1 b)(a _1 b)“ 1 
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or as 


, . T.-l T , ,-l 

a + (a ) c c ■» a-bd c 

and 

(a^a+c^c) (a ^b) = -c^(d-ca ^b) 

or as 

[ (a T )“ 1 c T +bd"' 1 ]c - 0 

and 

ab + cca b + cd-cca b = 0 

or as 

(c T d+a T b)d~ 1 c =0 (27) 

and 

(c T d+a T b) - 0 (28) 

Satisfaction of Eq. (21) is sufficient for Eq. (27) and both necessary 
and sufficient for Eq. (28), so the Proposition is proven, 

INTERPRETATIONS FOR STRUCTURAL DYNAMICS 
Eq. (22) is obviously valid when M is block diagonal, since then b=c=0 
and Eq. (21) is trivially satisfied. This result is consistent with Eq. (14). 
Eq. (22) is also valid when 

c = -a and d = b (29) 

and this result has an important physical interpretation for structural 
dynamics . 

In assessing conditions for the equivalence of Eqs. (17) and (18), we 
must recall that the columns of $ are the eigenvectors of Eqs. (9) and (10), 
and hence the "mode shapes" of independent structural vibrations for certain 
boundary conditions. The columns of the partitions $ and $ (see Eq. (14)) of 
$ represent respectively those modes which are preserved and those which are 
abandoned in truncation. Although several truncation criteria are useful, 
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when the elastic structure is an appendage on a rigid body over which one 
must maintain precise attitude control but not precise position control (as 
in the spacecraft problem) the essential criterion is the truncation of inodes 
which have no influence on the attitude of the primary rigid body to which 
the appendage is attached. If the mode shape is symmetric with respect to 
the primary body, as in Fig. 1, then clearly it must be removed in truncation, 
while antisymmetric modes, as in Fig. 2, must be preserved. In the special 
case in which all modes are either symmetric or antisymmetric, it becomes 
possible to formally justify removal of the former by the proposition of the 
previous section, as will be shown. 

Eq, (8) indicates that Eq. (13) must in greater detail be representable as 



i = n 

i j_ 

i ® . 


(30) 


where A and A are diagonal matrices of eigenvalues in the preserved and 
deleted categories, respectively, while <j) and cj> establish corresponding mode 
shapes for the position variables in q. Application of the criterion in 
Eq. (21) to Eq. (30) produces the requirement 

X? T <$ + - 0 (31) 

for the equivalence of Eqs, (17) and (18), and this is for all practical 
purposes never satisfied, since it reduces to the requirement 

X X, - 1 (a = 1, . . . ,N; $ = 1 6n-N) (32) 

However, if the problem is restructured so that Eq. (8) is replaced by 

Qi ... q 6n q 6n 3 T (33) 

in which q^,...,q an and <l 3 n+ ^» * * • »^6n 3re coor< ^ nates defining parallel 
motions of appendage nodes located symmetrically with respect to a central 
rigid body (as would be possible for the system in Figs. 1 and 2, but not 
generally so), then instead of Eq. (30) we can write 
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I = - 4 - 


where the left partitions (comprising $) contain the antisymmetric modes and 
the right partitions contain the symmetric modes. Now Eq. (29) is satisfied, 
and (with equations restructured as required by Eq. (33)) Eqs. (17) and (18) 
are identical, since inversion and truncation are commutative. 

SPECIAL CASE WITH D* - A* = 0 

When energy dissipation is ignored, and any rigid nodal bodies in the 
appendage model are in the nominal state spinning about principal axes of 

ft 

inertia, then D = A = 0 [1] in Eq. (1), which becomes 

Mq+Gq+Kq=L (35) 

In the undamped stable case of interest here all eigenvalues from Eqs. 
(9) and (10) are imaginary, and may be designated in conjugate pairs as 
±i0 a f° r a = l,...,6n. The corresponding eigenvectors are however complex, 
suggesting the designation 

+ ir a (36) 

where the form of Eq. (8) requires 


If in Eq. (35) we have L =0, we can without contradiction substitute 


q - 4>z + y a _1 z 


(38a) 


q = -yaz + tpz 


(38b) 


q = -yaz + \pz 

into the homogeneous part of Eq. (35) to obtain 


(38c) 


M (-yaz + ij/z) + G (-yaz + tpz) + K (ipz + yaz) = 0 
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( 39 ) 


t » t i i . i 

M i|/z + (-M ya + G ip + K y a )z + (-Gyo + K ip)z = 0 

A homogeneous solution of Eq. (35) is 

io t 

q h = (<J> + iY )e (40) 

which when substituted yields 

m'^ 01 + iy a ) (-°q) + c’of + iy a )(i a a ) + K(^ a + iy a ) = 0 
or 

-M i|/*a a 2 - G y a a a + K ip a = 0 (41a) 

and 

-m’y 0 ^ 2 + g'^o + k'y 01 = 0 (41b) 

Ct (X 

and Eqs . (4 la) and (41b) considered for a = l,...,6n imply 
i 2 * * 

-M fa - G yo + K l() a o (42a) 

-m'yo + g\ + k'y 0~ X - 0 (42b) 

Thus with Eq. (42b) Eq. (39) reduces to 

i i 

M ipz + (-Gyo + K ip)z = 0 

T 

which after premultiplication by ^ becomes 

^ T M \\>z + (i]j T K \j) - \ji T Gya)z = 0 
X 

and with times Eq e (42a) this is 

^ T M*4»(z + cr 2 z) = 0 
T 1 

or with nonsingular 4> M ijj, 

z + a 2 z = 0 (43) 

This same result arises if Eq. (38b) is equated to the derivative of Eq. 

(38a), so that contradiction is avoided. The transformations in Eq. (38) 
are not acceptable in the inhomogeneous case of Eq. (35) of primary interest 
here, however, because these transformations directly imply Eq. (43). Instead 
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Eq. (35) must In the inhomogeneous case be written in first order form, as 
in Eqs. (9) and (10). In the latter case, one can accomplish the desired 
transformation to uncoupled (pairs of) scalar equations with the real 
transformation 

Q - PZ (44) 


where 


P A 


—t 

1 X 

-yo~ 

l W 


(45) 


,-l 


and P exists if the eigenvectors of the system are independent, as required 

for stability in this undamped case. Substituting Eq. (45) into Eq. (10) 

* * _1 

with D = A =0 and premultiplying by P produces 


Z = 


„ 2 _ 

*-o" 


a 

'O' 


Z + P -1 L 


(46) 


as may be confirmed by considering the homogeneous case and defining 


Z = 


-1* 

0 Z 


Eq. (46) offers a substantial advantage over Eqs. (12) and (16) , which 
involve complex numbers in Y, A, and $ , although Eq. (12) possesses an 
advantage in requiring the inversion of only a diagonal matrix. The 
inversion of P in Eq. (46) becomes particularly critical when (as required 
for practical structural dynamics) this inversion is to be followed by, or 
precede, coordinate truncation. We must once again face the critical question 
raised in this paper, "Are truncation and inversion commutative?” More 
explicitly, we must determine whether or not the truncation of Eq. (46) to 


Z = 


0 

— 1 

a 

-o’ 

L, 

0 


Z + (P -1 )L 


(47) 


provides the same result as would be obtained by substituting 
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z 


(48) 


= PZ A 


Ljl 


. -Y o 

! P . 


t » 

into Eq. (10) with D = A = 0 and premultiplying by 


P fc A (P T P) -1 P T 


(49) 


This issue was not confronted in [2], where an affirmative answer was pre- 
sumed. This presumption requires for its validity that 


and 




(50) 


(51) 


Eq. (50) is proven in an appendix of the referenced report which underlies 
[1], hut the validity of Eq. (51) has not been examined previously. With 
the proposition presented as the central result of this paper (see Eqs. (21) 
and (22)), it becomes evident that Eq. (51) is correct if any only if 


-a T Y T i^ + = 0 (52) 

and this is an untenable hypothesis, Eq. (51) must be rejected, and the 
very appealing real transformation represented by Eq. (48) cannot be accepted 
as the formal equivalent of the truncation of the rigorously valid Eq. (46). 


CONCLUSIONS 


The necessary and sufficient conditions presented in Eq. (21) for the 
commutativity of coordinate truncation and transformation matrix inversion 
(see Eq. (22)) are quite severe. Two of the three transformation procedures 
considered here for flexible appendages on constantly rotating rigid bodies 
are essentially destroyed by the noncommutativity of these operations, since 
for systems of realistic dimension it is not feasible to perform the necessary 
matrix inversions prior to truncation, as an honest interpretation of the 
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mathematics of the problem seems to require. This discouraging result is not 
yet a definitive conclusion, because truncation is at best a process of ap- 
proximation. It may yet be true that transformations such as Eqs. (20) and 
(48) give satisfactory results in many cases; at this point however there 
appears to be no acceptable evidence supporting this possibility. 

The only one of the three transformation procedures considered here 
which survives the commutativity test is that culminating in Eq. (15). This 
procedure involves no inversion except for a diagonal matrix, but it is a 
transformation to complex coordinates, and it generally requires not only 
the selected system eigenvectors in $ but also the corresponding adjoint 

__i « » 

eigenvectors in $ . In the important special case for which D = A =0 
this procedure is much simplified by the relationship [1] 

?’ = tf (53) 

f 

but unless the appendage is nonrotating (so that G =0) the transformations 
involve complex numbers. 
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CHAPTER III 


MATCHED ASYMPTOTIC EXPANSION 
MODAL ANALYSIS OF ROTATING BEAMS* 

Andre D. Collnt 
and 

Peter W. Likinstt 
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ABSTRACT 

Modal analyses are presented, for the transverse vibrations of a uniform 
Euler-Bernoulli beam, rotating about an axis orthogonal to the beam. Results 
are obtained for a cantilever beam emanating radially from a base which is 
rotating but not translating in inertial space, and for a beam with built-in 
ends spanning the diameter of a ring which is rotating about its inertially 
fixed symmetry axis. The equations are formulated in terms of a small 
parameter e which is proportional to beam stiffness and inversely related to 
mass density, length, rotation rate, and the pretension of the diametral 
beam. The resulting singular perturbation problem is analyzed by the method 
of matched asymptotic expansions, employing asymptotic expansions for the 
central region and for boundary layers at either end. Results are expressed 
literally, with sample numerical examples. 
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MATCHED ASYMPTOTIC EXPANSION 
MODAL ANALYSIS OF ROTATING BEAMS 

Introduction 

The rotating elastic beam has technological significance in the develop- 
ment of helicopters and spin-stabilized space vehicles, and it has been analyzed 
extensively. Because the linearized partial differential equations of 
vibration have coefficients depending on the independent spatial variable, 
closed form exact solutions are not available, and recourse to approximate 
methods is necessary. Most investigators have adopted numerical analysis 
procedures, obtaining results that are limited to the particular parameter 
values considered. Yntema 1 provides a good illustration of the extensive 
application of Galerkin's method to many specific cases, expanding the 

modes sought in terms of the nonrotating beam modes; Renard and Rakowski 2 

3 

and Hughes and Fung use similar numerical methods to address specific 
problems of interest in their work. 

In this paper we seek literal expressions for approximate representations 

of natural frequencies and mode shapes, using the method of matched asymptotic 

4 

expansions. Modal analyses of rotating structures have been developed 
previously by Boyce and Handelman 5 and by Abel and Kerr . 6 In the former 
paper the authors generate zeroth order solutions for a rotating beam with a 
tip mass, applying a method developed by Moser 7 in which the linearly 
independent solutions are taken in the form of B(x,ry) expln - 1 h(x) ] , where r| 
is of the order of magnitude of the small perturbation parameter of the problem 
and B(x,q) and h(x) are functions to be determined. In the paper by Abel and 

g 

Kerr the method of asymptotic expansions is applied to a rotating cable- 
counterweighted space station for cables with small flexural rigidity. Match- 
ing of central and boundary layer expansions is employed to establish those 
constants of integration not established by boundary conditions. 

, ,63 Preceding page blank 
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In what follows the rotating Euler-Bernoulli beam vibration equations are 
presented and a central solution expansion developed, and then boundary layer 
expansions are developed and matched for the boundary conditions appropriate 
for rotating beams of both "radial" and "diametral' 1 configuration. 

Rotating Beam Equations 

The small transverse vibration W(£,t) of an Euler-Bernoulli beam with 

constant flexural stiffness El and constant mass per unit length ]i, under 

8 

the external axial load per unit length P(£), is characterized by 



In what follows we address two distinct problems of beam vibration: the 

radial rotating beam (Fig. 1), and the diametral rotating beam (Fig. 2). In 
each case the beam is attached to a body B which has a prescribed constant 
inertial angular velocity Qe^ with the £,n,5 cartesian coordinate system axes 
fixed in B and its origin. 0 fixed both in B and in inertial space. In each 
case the beam supports are built into B, constrained against translation and 
rotation relative to B. Thus the boundary conditions are as follows: 


2 3 

Radial Beam: W(0,t) ® “ (0,t) = (L,t) = (L,t) » 0 

Diametral Beam: W(L,t) = W(-L,t) * ||- (L,t) = — (-L,t) = 0 


( 2 ) 

(3) 


In each case the rotation induces axial strains in the beam, with 
accompanying change in transverse stiffness characteristics . As has been 

g 

formally established by the use of nonlinear strain-displacement equations, 
one can for the small strain case continue to represent small transverse 
vibrations by Equation (1), incorporating into P(£) the "centrifugal force," 
Thus for the radial beam 
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For the diametral beam, a pretensile load T is required to avoid instability 
of the solution W(£,t) = 0; thus the effective axial "force” is 


2(0 « T 



VtthdZ - T - yfl 2 £ 2 /2 


(5) 


with 

T > yfi 2 L 2 /2 (6) 

2 

By defining k such, that 
2 2 2 

T A k pDT/2 

one can write Equations (4) and (5) together as 

PCO - un 2 (k 2 L 2 -£ 2 )/2 (8) 

2 2 
where for the radial beam k - 1, and for the diametral beam k > 1. 

Equation (1) can now be written in a form applicable to both radial and 

diametral beams as 


El 




2 


3_ 

9? 


[(k 2 L 2 -? 2 ) 


3W 

3? 


*1 & 2 W 

J dt l 


Representation in terms of dimensionless quantities is accomplished by 
introducing 

w A W/L 
x A £/ (Lk) 

T A Jit 

2 

and dividing Equation (1) by yfi L/2, to obtain 


e 


3 4 w 

9x 4 


8 


[a-* 2 ) U}+ 2 


2 



(9) 


CIO) 


where 



Equation (10) admits the separable product solution 

w(x,t) * <f)(x) a(t) (12) 

with 

e<f> IV - [(l-x 2 )<j> , j’ - \ 2 <p = 0 (13) 

and 

a + ~ \ 2 o - 0 (14) 

2 

for some constant A , with prime denoting d/dx and dot denoting d/dT. Those 

2 

discrete constant values of A for which Equation (13) has nonzero solutions 
are the system eigenvalues . 

In terms of the separated dimensionless quantities in Equation (13), 
the boundary conditions become for the radial beam (for which k = 1) 

<K0) = 4> f (0) - 4> T? C L) - 4> m (1) - 0 (15) 

from Equation (2), and for the diametral beam 

<)>(k 55 <K-k ■ $ ! (k - <t>’ (-k ^) « 0 (16) 

from Equation (3). 

The objective of this paper is to obtain explicit solutions for the 
modal functions <f>(x) in Equation (13) as asymptotic expansions in e, for 
boundary conditions as in Equations (15) and (16) * This process must also 
yield asymptotic expansions for the system eigenvalues, which produce natural 
frequencies of vibration from Equation (14) and the time normalization 
t « fit as 

to * ( 17 ) 

The solution technique employed is the method of matched asymptotic expansions, 
as presented by Cole. Because in the limit as e 0 Equation (13) becomes 
only a second order differential equation, the solution of which cannot 



satisfy the four boundary conditions in Equation (15) or Equation (16), this 
problem is classified as singular . The solution will be obtained as the 
combination of an asymptotic expansion valid in the central region of the 
beam, obtained from the reduced order differential equation arising in the 
limit e + 0, and asymptotic expansions valid near the boundaries of the 
beam. 

Central Solution Expansion 

For the central region of the beam, we represent the dependence of the 
solution of Equation (13) upon the small parameter e by the asymptotic 
expansion 


<Kx,e) = h Q (x) + v 1 (e)h 1 (x) + v 2 (e)h 2 (x) + ... 
where v^e) is for i = 1,2,... an asymptotic sequence to be established. 


(18) 


Similarly the eigenvalues represented by A in Equation (13) are expanded as 
A 2 = Aq + K 1 (e)A 2 + < 2 (e)A 2 + ... (19) 


for some asymptotic sequence K^(e), i = 1,2,... 

Substitution of Equations (18) and (19) into Equation (13) and considera- 
tion of the limit e -*■ 0 requires 


I(l-x 2 )h^' + A 2 h 0 = 0 

Equation (20) is a Legendre equation of order n, where n(n-KL) = The 


( 20 ) 


general solution of Equation (20) is given by 


10 


where 


h A = AU (x) + BV (x) 
0 n n 


. ‘ . V . , A 0 2 VV 6 > 4 
V X) - 1 ' 2 \* + — 4! x " 


(21) 


( 22 ) 


and 


al-2) , (A^-2)(a1i2) 

■'.« A * - -IT- * + -V 


x - ... 


(23) 
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2 

The Legendre functions U n (x) and V n (x) converge for arbitrary A Q over the 

range -1 < x < 1, and either ^(x) or ^(x) terminates for those discrete 
2 

values of Ag for which 

Aq = n(n+l) n = 0,1,2,... (24) 

with U (x) terminating for n even and V (x) terminating for n odd. The 
n n 

terminated series can then be expressed in terms of a Legendre polynomial 
P n (x) by 

U (x) - P (x)U (1) n = 0,2,4,... (25a) 

n n n 

or 

V (x) * P (x)V (1) n = 1,3,5,... (25b) 

n n n 

and the non terminating series can be written in terms of Legendre functions 
of the second kind, Q^(x), which are unbounded in the limit x = 1. These 
observations are particularly significant for the radial beam, as indicated 
by the boundary conditions in Equations (15) and (16). 

In addition to the zeroth order approximation given by Equation (20) , 
one obtains higher order term equations from the substitution of Equations 
(18) and (19) into Equation (13). If no restrictions on V^(e) or ^(e) 
were stipulated, the first equation involving terms above the zeroth order 
would be 

eh 0 IV - Vj.CeMO— x 2 )h^] ' - v^eM^Ce) = ^(eJA^ (26) 

Further progress requires knowledge of the order relationship among K^(e), 

•j^ 

V^(e) and e. The possibility that K^(e) = °(v^(e)) and K^(e) = o(e) is 

2 

rejected, because Equation (26) then indicates that A^ = 0. Similarly the 
possibility that £ = o(v^(e)) and E = o(k^(£)) is rejected because of the 
implication that hg = 0. Subject to later confirmation, we assume that 
v^(e) = o (e) and ic^(e) = o(e). Thus Equation (26) becomes either 

By notational convention a( ) = o(b( )) implies ***■“ (—) = 0. 
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(27a) 


[(l-x 2 )!^] + A 2 ^ = -A 2 h Q 


I^i (e) - 


or 


t Cl-x 2 )hj] ’ + A 2 h x 


0 [v 1 (e) = oC^Ce))] (27b) 

In order to proceed further with this calculation we require v^(e), which 
must be established by the requirements of matching the expansions which are 
valid at the boundaries. These expansions are different for the radial beam 
and the diametral beam, and these cases must be treated separately. 

Radial Beam Inner Expansion 

In the neighborhood of the boundary x = 0, the term in Equation (13) 

becomes important, and the previous solution is invalid. In order to develop 
a solution which is valid in the "boundary layer" near x = 0, we introduce a 
stretched coordinate 

x A x a ^(e) ( 28 ) 

assuming o(e) ■+■ 0 when e -*■ 0. Into the original Equation (13) we now substi- 
tute the asymptotic expansion 

4>(x,e) = p 0 (e)g 0 (x) + P 1 (e)8j t (x) + P 2 (£)g 2 (x) + ••• (29) 

where p^(e) constitute an asymptotic sequence for i 0,1,2,..., to obtain 
( (e) d 4 g 0 (x) P 1 (e) d^g^x) 


P 0 

L 4 


(e) dx 


o 4 (e) 


dx 4 


+ . . . 


- (1-a (e)x ) 


2 | i-q^) d 2 g Q (x) p-^e) d 2 g 1 (x) 


fo 

.a 2 


(e) dx 


a 2 (e) dx 2 


■] 


+ 2 a(e) 


t p Q (e) dg Q p 1 (e).d 8 1 (x> . ■ 1 

V(e) dx" + o(e) ir “ + 


- Aq [Pq(e) g Q (x) + p-^e) g x (x) + . ..] 


- fc 1 (e)A^[p 0 (e) g Q (x) + P 1 (e) g L (x) + ...] = 0 
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The coordinate scaling factor o(e) is chosen so as to maintain the 
highest order derivative and an additional term or terms in the lowest order 


approximation of this expansion. The choice 



then produces the dominant boundary layer equation 

d 4 g Q (x) d 2 g Q (x) 

— Zk " 0 

dx dK 

which after two integrations provides 


(30) 


(31) 


d 2 So® . -i , . i 

— - c o e + c i e 

dx 


(32) 


The inner expansion in Equation (29) in the limit as i + “ must match the 
central expansion in Equation (18) in the limit as x ■> 0, so that in 
Equation (32) must be zero. Two more integrations then provide gg(x), which 
with the boundary conditions 

d 8nl 


&o (5) 


= 0 


(33) 


x=0 


dx 


x=0 


must be 

g Q (x) = C Q (x - 1 + e -X ) (34) 

In order to match the approximation of the inner expansion (Equation (29)) 
provided by Equation (34) with the approximation of the central expansion 
(Equation (18)) provided by Equation (21) we can introduce an intermediate 
limit defined by the coordinate measure 


x n A xr 1 _1 ( £ ) 

1/2 

with n(e) ■* 0 as e -*■ 0 and t = o(n(e)), so that in the limit as e -*■ 0 with 

~ -1 /2 

x ^ fixed, x = n(c)x^ •+■ 0 and x = ri (e) e x^ -*■ 00 . Matching then requires 
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Lim 

{h (nx ) + V- (e)h. Cnx ) + 

x fixed u n 1 1 n 

n 




+ . . . } 


( 35 ) 


Expanding hg(x) and h^(x) as Taylor series about the origin x = 0 and substi- 
tuting Equation (34) for gg(x), we find 

Lim , 

e-^0 (h ft (0) + nxh' (0) + ... + V 1 (e)h. (0) + V (e)nx h.(0) + ... 

x u nu xx xrii 

n fixed -1/2 

I/O ^ 

- p 0 (e)C Q (ne ' x n -1 + e n ) +...} = 0 (36) 


Equation (36) is satisfied by 


P 0 Ce) 

(e) 


= e 1/2 


-e 1 ' 2 


v°> 

h 0 (0) 

h a (0) 


0 

C 0 

-c 


0 


(37) 

(38) 

(39) 

(40) 

(41) 


which conforms to the assumption v^Ce) = 0 (e) preceding Equation (27a). 

With Equations (39) and (40) we can return to Equations (21) - (23) to 
conclude that 


A = 0 and B = (42) 

Thus first approximations are fully established for both inner expansion 
(Equation (29)) and central expansion (Equation (18)), in terms of the constant 
C Q . Moreover, we can now recognize that the first approximation of the 

eigenvalue expansion in Equation (19) must belong to the set of discrete 

\ 

values for which the Legendre function V n (x) in Equation (23) is a finite 
series; otherwise V n (l) and 0(1) become infinite. Thus Equation (24) must 
be satisfied for n odd, so that 
2 

Aq = 2, 12, 30, ... n(n+l) for n = 1, 3, ... ( 43 ) 
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and the natural frequencies identified in Equation (17) are approximately 

u) 1Q = ft, ftv^6, ftyTif ^[n(n+l)/2] 1 ^ 2 for n = 1, 3, ... (44) 

The next step is the determination of h^(x) in Equation (18). We make 
use here of the orthonormality of the two different solutions 4> a (x,e) and 
ij>g(x,£) corresponding to the distinct eigenvalues X^ and Xg, which requires 



d> cf> _d m 6 „ 
y a Y B x «B 


(45) 


where 6 0 is the Kronecker delta. Equation (45) can be proven for a f 3 by 

Otp 

substituting <f> for 4> in Equation (13), multiplying the equation by and 

C6 P 

integrating twice by parts, using Equation (15) to obtain 


' c *a+e 


1 

+ E 
0 "0 


1 + / a- x 2 )<t’^dx - X 2 / * a * e - 0 


and then repeating the process with a and B exchanged, finally subtracting 

2 2 

these results from one another and noting that X^ ~ / 0* The validity 

of Equation (45) for 6 = 0i is freely prescribed, since there is a free 
constant factor to be prescribed in any solution to Equation (13) . 

If <J) and 4>„ are asymptotic expansions in e, then Equation (45) equates 

Ot p 

an asymptotic expansion to 6^. For this equation to be valid for any e, 

the zeroth order term must be 5 and all other terms must be zero. We 

a3 

have established that in the inner boundary layer the first term is of order 

1/2 

e (see Equations (29) and (37)), and moreover the boundary layer "thickness 1 ' 

1/2 

is of order e (see Equation (30)). Thus the contribution of the inner 

expansion to Equation (45) is limited to terms of order e and above. If we 

now assume (subject to later confirmation) that the outer expansion contri- 

1/2 

but ion to Equation (45) is also of lower order (higher power) than e , then 
we can substitute the central expansion (Equation (18)) alone into Equation (45) 
to obtain, for a single mode (3 ■ a) 
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h 2 (x)dx = 1 


^ h 2 (x)dx = 1 
£ h^x^^dx = 0 


( 46 ) 


(47) 


The combination of Equations (21), (42), and (46) yields 


' [/ lv n (x>te ] 


-1 


(48) 


Equation (47) is useful in resolving the choice between Equations (27a) 
and (27b); in the latter case h^(x) is simply a multiple of hg(x), and 
Equation (47) is impossible. Thus we now conclude (noting Equation (38)) that 


«1 (e) = v x (e) = e 1 ^ 2 

and h x (x) must satisfy Equation (27a) and Equation (41). 

2 

Equation (27a) involves A^, which can be determined by multiplying 
Equation (27a) by h^(x) and integrating by parts from 0 to 1, to obtain 

,1 


(49) 


(l-x 2 )hjhgdx 


a -* 2 ) h ’ h 0 | o - y 1 

■o /w* '-•it h o 


+ A 


[dx 


0 *** t) 

or, with Equations (39) * (46), and (47) and another integration by parts. 


- (1-x )h Q h 1 


■ * f 


[(l-x 2 )hQ}’h 1 dx - - A 2 


[0 *'0 

2 

Substituting - A^h^ from Equation (20) into the final integrand and utilizing 
Equation (47) again, we find 

A 1 “ " V 0 ) h i (0) 


With Equations (40), (41), and (48), this becomes 


73 



(50) 


A 1 ‘ [/ V n (K)dx ] 


More specifically, corresponding to the first two inodes defined by Equation (43), 

(51a) 


A ^- 


/ V- (x)dx = / x dx 

L-6 1 J L o J 


- 3 


and 


, ■ [/’ •;■»■] “ 

[f 


(x - X 3 ) 2 dx 


r- 


63 


= 15.75 


(51b) 


2 2 

(It may be noted that Equations (43) and (49) for Aq and could alterna- 
tively have been obtained from Equation (45).) 

2 

Finally we can return to Equation (27a) with A^ from Equation (49) and 
seek h^(x) from that combination. The homogeneous solution is of course 
merely some factor k^ times h^, and since h^ (0) = 0 this leaves the particu- 
lar solution h^(x) to satisfy the boundary condition on 1^(0) in Equation (41), 
Thus 


h^(x) = h^(x) + kjhg 


(52) 


where from Equation (47) 

0 = J' h 1 (x)h Q (x)dx * J h 1 (x)h {) (x)dx + J' h^dx 


so that 


- I h l (l!)h 0 


(x) dx 


(53) 
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We can proceed without a more explicit expression for h^(x) to accomplish 
the necessary matching with an outer asymptotic expansion to be obtained for 
a boundary layer near the free end of the beam. 

Radial Beam Outer Expansion 

At the free end of the radial beam in Figure 1 we must again satisfy 
boundary conditions which cannot be reconciled with the central expansion. 

We therefore introduce another boundary layer coordinate 

x* A (l-x)a -1 (e) (54) 

with a(e) 0 when e ■+ 0. 

Into the original Equation (13) we now substitute the asymptotic expansion 
<K*.e) = 6 0 (e)f 0 (x*) + 6 1 (e)f 1 (x*) + ... (55) 


where 6 . (e) constitute an asymptotic sequence for i = 0, 1, 2, .... to obtain 

I V £ > d V* 

e 


r 6 

i 

.a 


Q (e) d 4 f.(x ) 6. (e) d 4 f.(x*) 

+ -4 =— — + 


4 ft) dx* 4 


a 4 (e) dx* 4 


• i • 


- [l - (l“ 0 t(E)x*) 2 


][) 


(e) d 2 f Q (x*) 6 x (e) d 2 f 1 (x*) 
— + — — + 


(s) dx 

* 


*2 


2 , . *2 

a (e) dx 


- 2 


* r<5 0 (e) df Q (x ) 6. (e) df (x ) 1 

a-c, (e)x > -0-^-0 + ^1_ + ... 

L dx dx 


- [6 Q (e)f 0 (x ) + ) + ...] 


] 


- ^(e) A 1 [6 0 (e)f 0 (x ) + 6 1 (e)f 1 (x )+...]= 0 

The choice 
ct(e) = 

produces from Equation (56) the dominant boundary layer equation 

d 4 f Q (x*) * d 2 f (x*) df (x*) 

- 2x — 2 — “r* = 0 


dx 


*4 


dx 


*2 


(56) 


(57) 


dx 


(58) 



Boundary conditions at x * 1 must be satisfied; from Equations (15) and (54) 

“ 0 (59) 


3 *. 

d f Q (x ) 


dx 


*3 


* 

x =0 


and 


2 * 
d f Q (x ) 


dx 


*2 


A 

x =0 


* 0 


(60) 


Integrating Equation (58) and enforcing Equation (59) produces 


3 * 

d f Q (x ) 


dx 


*3 


* df 0 (x > 
< — * — 
dx 


(61) 


which can be recognized as an Airy equation in the variable df^(x )/dx , with 
a boundary condition established by Equation (60), and solutions available'*'^' 

{te ft 

as a linear combination of the Airy integrals Ai(x ) and Bi(x ). However, 

Vc Jc * 

Bi(x ) cannot appear in the solution for dfg(x )/dx because it grows expo- 

A 

nentially with some power of x , and matching with the central solution 

ft ft 

expansion in the limit as x -*■ » would be impossible. Moreover, Ai(x ) 

cannot appear In the solution because it cannot meet the boundary condition 

at x * 0 imposed by Equation (60) . Thus the required solution to Equation 

(61) must be zero, implying 


f o (x * ) 


(62) 


a constant . 

In order to match the outer expansion in Equation (55) with the central 
expansion in Equation (18), we again introduce an intermediate limit defined by 
the coordinate measure 

= (1 - x ) n ( e )"* 1 

1/3 

with n(s) + 0 as £ ■+ 0 and e « oCn(e)), so that in the limit as e + 0 with 
x n fixed, (1-x) = x n n(e) + 0 and x* = (l-xje" 1 ^ 3 * x^rKeH -1 ^ 3 «. Matching 
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in this limit then requires (in view of Equation (38)) 


Lim 


e+0 {h n (l-nx ) + (1-nx ) + ... 

x fixed u n ± n 

T) 

- + 6 1 (e)f 1 (e~ 1/3 n x ri ) + ...}= 0 

Expanding h^, and h 2 as Taylor series about x^ ■ 0 yields, with Equation (62), 

t 1 tl ry J 1 m -y •, 

X £ f ixed {h ° (1> ' h o (1)r|X n + 27V 1)n x n + 3i\ (1)n x n + 

n 

+ e 1/2 h 1 (l) - E 1/2 h^(l)nx n + ... +v 2 (e)h 2 (l) + ... 

- 6 0 (e)D 0 - 6 1 (e)f 1 (e _1/3 nx Ti ) - 6 2 (e)f 2 ( G “ l/3 nx n ) + ...}-= 0 (63) 


Matching requirements suggest the tentative adoption of the following equations: 
v 2 (e) = e (64) 


VD - d 0 


6 0 (e) - c - 1 


^(e) = e 


1/3 


<$ 2 (e) ■ £ 

<$ 3 (e) = e 
6 4 (e) - e 


1/2 


2/3 


5/6 


<$ 5 (e) = e 


(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 


and so forth. These values are compatible with the assumption made prior to 
Equation (46), and provide the means to return to Equation (56) to establish 
differential equations to be solved for f^(x ), f 2 (x ), .... With Equations 
(49), (57), and (66) - (71), Equation (56) can be rewritten as follows: 
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-1/3 ^fo- 


lds* 4 

fA 

L** 4 ' 

+ £ 1/6 f 


dx 


*2 


df. 


dx 




+ e 


* ^ ^1 *2 ^ ^0 ^1 * 

2x —£ + x -jg - 2 -± + 2x 

dx dx dx dx 


- A 2 f I 

« * A o 0 

dx J 


d 4 f. 


2 * ^ 

LT?4- 2X ~ 7~*2 ~ 


df 


dx 


dx 


] 


+ £ 


1/3 


+ e 


1/2 


d 4 f„ 


* d2f 3 


*2 d2f l 


. *4 

dx 

dx* 2 

*2 

dx 

d ' f 4 

A 

O.* *+ 1 

*2 d2f 2 


df 3 * «! 2 

-2 A + 2x * \)^1 

dx dx u 1 


] 


*4 


dx 


*2 


dx 


*2 


df 4 * df 2 a 2 

2 * + 2x * ^o f 2 


+ e 


2/3 


dx 

df 


dx 


- A ? f o] 


r^5 * ^ ^5 2 ^ ^3 ^^5 * ^3 2 

[”*4 ‘ 2x ~* 2 +x T "* 2 ' 2 T "* + 2x 7 * ‘ A 0 f 3 
L dx dx dx dx dx 


] 


+ ... = 0 (72) 

2/3 

assuming that < 2 (e) = o(e ) in Equation (19). 

Tha validity of Equation (72) for any e requires that each of the ex- 
pressions in square brackets be zero. The first such equation has already 
been recorded as Equation (58), and its solution for the required boundary 
conditions has been given as Equation (62). The second bracketed expression 
in Equation (72) then simplifies to 


d \ * d2f l df l 2 

dx* 4 ' 2X dx* 2 ' 2 dx* ’ V0 < 73 > 

For f^tx ) and terms of higher index the boundary conditions on the free end are 


d 2 f . 


dx 


*2 


A 

x =0 


- 0 


j ~ 1 » 2, 3, ... 


(74) 


d 3 f . 
J. 


dx 


*3 


A 

X =0 


0 


j = 1, 2, 3, .. . 


(75) 
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Equation (73) has the same homogeneous solution as Equation (58), and in 
addition the particular solution 
— , * 1 2 * 

f l (x ) “ " J A 0 D 0 X (76) 


Satisfaction of the boundary conditions again eliminates the Airy integrals 
in the general solution, so that Equation (76) with Equation (65) provides the 
total solution. 

Now we can evaluate the validity of Equation (67) by returning to 

\ i 

Equation (63) to examine the matching of the f^ term and the h^(l) term, 
which together become 

-h'(l)nx n - 6 1 (e)f 1 (e“ 1/3 nx Ti ) - 

[ -VD + | e 1 / 3 A^h 0 (l)e _;L/ 3 ]nx n - [-h'(l) + \ A^U)]^ 

The expression in brackets is indeed zero, as required by Equation (20) 
evaluated at x = 1, and these terms are properly matched. 

The third bracketed expression in Equation (72) has the same structure as 
the first, and admits a solution reduced by boundary conditions to a constant, 
which for matching in Equation (63) must be h^(l); thus 

f 2 (x*) = 1^(1) (78) 


The bracketed term multiplied by e 
Equations (76) and (65) to provide 


1/3 

in Equation (72) combines with 


* d2f 3 

- 2 X — d - 2 


df. 


dx 


*4 


dx 


*2 


dx 


(2 - Aq) Ajh 0 (l)x*/2 


( 79 ) 


with the particular solution 



and the general solution involving Airy integrals Ai and Bi as previously. 
Again Bi must be rejected for its exponential growth, which precludes 
matching for x but here Ai cannot be rejected, and the general solution 

to Equation (79) becomes 

* 

f 3 (x*) = A 2 (A 2 - 2)h 0 (l)x* 2 /16 + D 1 jf AiCOd? + D 2 (81) 

Equation (74) imposes on solution (81) a boundary condition providing 
A 0< A 0- 

1 8Ai'(0) (82) 

' 11 
with Ai (0) from tables. 

II 

Matching in Equation (63) requires that the term involving h Q (l) sum to 

-1/3 

zero in the limit with the unwritten term - 6 3 (e)f.j(£ nx^) . That portion 
* 

of f 3 (x ) involving and decays exponentially, so that matching requires 
only 

1, "... N 2 2 2/3. 2, .2 ... ... -2/3 2 2... 

jjhgCOn x^ - e A o^O " 2 ^ h oU) £ ^ x n / 16 

- [8h”(l) - (A 2 - 2)A 2 h 0 (l)]n 2 x n 2 /l6 = 0 

This equality holds since the expression in brackets is zero by virtue of 
Equation (20), which may be written in terms of Taylor series about x m 1 as 

[- 2 h'(l) + A 2 h 0 (l)] + [4hQ<l) + 2 h’(l) - A 2 h' ( 1)](1 - x) 

+ [-3hQ (1) - 3hg(l) + y A 2 h 0 <l)Kl - x ) 2 + ... = 0 (83) 

so that each of the bracketed expressions is zero. 
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Solving the three equations implied by Equation (83) simultaneously 
produces 

V 15 -J A o h o (1) (84 > 

h > - J <A 2 -2)A 2 h 0 <l) (85) 

h o a) -h < A o- 6 > «o- 2 )A o h o (1 > (U) 

The next step is the solution of the equation involving the coefficient 

1/2 

of e in Equation (72). With Equations (78), (62) and (65), this equation 
becomes 


d4f 4 * d2f 4 

- 2x — rf - 2 


df 


*4 


*4 


J = A^Cl) + A 2 h rt (l) 


0 0 


(87) 


dx dx ~ r dx 

The boundary conditions in Equations (74) and (75) limit the solution of 
Equation (87) to 

f 4 (x*) - - \ (AQh 1 (l) + A^h o a)) X * 


( 88 ) 


X 1 

Matching of the f^(x ) term with the h^(l) term in Equation (63) requires 
- e 1/2 h’(l)nx n + ~ e 5/6 (A 2 h 1 (l) + A 2 h 0 (l))e' 1/3 nx n 
« t-h’ (1) + \ (A 2 h x (l) + A 2 h 0 (l))]e 1/2 nx Ti - 0 


which is assured by Equation (27a) with x “ 1. 

To complete the projected treatment of the outer boundary layer, we set the 
last bracketed expression in Equation (72) to zero, noting Equation (81). 

Since our concern here is not with an exact expression for f^(x ) but only with 

jlf sV 

enough information about f^(x ) as x ■* 00 to establish the matching and verify 

& 

Equation (71), we write the equation for f^(x ) as 


d 4 f , d 2 £ df 

-s|- 2 * — *f-- 2 — i 

dx dx dx 


(Aq- 6)(A 2 -2) h 0 (l)x* 2 + TST 


(89) 
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where TST represents transcendentally small terms which have no influence on 
the matching process. Equation (89) is easily integrated, to obtain a result 
presented here as 


df 5 -(A 2 -6)(A 2 -2)A 2 h 0 (l)x* 2 1 d 3 f 5 

= 96 + “* *3 

dx 2x dx 


+ TST 


2x 


Two differentiations yield the approximation 


(90) 


d 3 f r 


dx 


*3 


<*0-M«o- 2H oV 1 > *-3, 

48 + 0(x > 


which in turn permits Equation (90) to be written as 

df 5 (A 2 -6)(A 2 -2)A 2 h 0 (l)x* 2 


dx 


D 3 + 


96 


(A 2 -6)(A 2 -2)A 2 h 0 (l) 
48 “ 


] 


2x 


„ *-4 . 

+ 0(x ) + TST 

Integration now yields the approximation 


, _ (A 2 -6)(A 2 -2)Ajh 0 (l)x* 3 

V* > 288 


L D 3 + 


(A 2 -6)(A 2 -2)A 2 h 

-43 


— ] + D i 


+ 0(x* 3 ) + TST 

The constants and must for matching be given by 

3 ' " 48 


(91) 


^ O jL ij 

The notation X = 0(x ) implies that Lim Xx 

* 

x ■* » 


is bounded 
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and in view of Equation (64) , 


° 4 ' V« 


Now matching in Equation (63) is assured by the identity 

37 h o (1>r|3:<3 n - ef 5 (£' 1/3 n Xl1 ) - 

K <1) , (A O- 6KA 0- 2)A 0 h 0a> 1 3 3 . 

I n x =0 


which follows from Equation (86) • 

Thus the outer expansion in Equation (55) has been completed to order e. 
Radial Beam Solution Summary 

In brief summary, the rotating radial beam in Figure 1, as characterized 

by transverse vibration equations (Equation (10)) with cantilever boundary 

conditions (Equation (15)) has natural frequencies given by Equation (17), 

2 

with X from Equation (19) and * (49) available as 


■ 4 ♦ c 1/2 4 


2 2 

where A Q is given for the sequence of inodes by Equation (43) and is given 
for all modes by Equation (50) and for the first two inodes explicitly by 
Equation (51). Although these results could be compared for specific spin 
rates to the computer-generated frequencies presented by Yntema and others, 
it is more revealing to compare to the formula for the first mode used by 
Vigneron: 

U>1 = w 2 R + 1.193A 2 

2 

where is the lowest natural frequency of the nonrotating beam, which is, 
in view of Equation (11), given by 

<4 = (3.515) 2 EI/(yL 4 ) = (3.515)Ve/2 
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Thus, in terms of our e, Vigneron uses 

w 2 = a 2 (1.193 + 3.515 2 e/2) (92) 

For comparison, our approximation is 

ojJ = Q 2 (l + 2.12i/e72) (93) 

These two approximations are protrayed graphically in Figure 3. More precise 
estimates indicate that Equation (93) gives better results for sufficiently 
small e (and yields the correct limiting case of the cable frequency when 
£ ■+ 0), but Equation (92) is more accurate for sufficiently large e. 

The mode shapes have been developed in this paper in terms of asymptotic 
expansions as follows: 

Inner boundary layer expansion 

4>(x,£) - C 0 £ 1/2 (£ _12 x - 1 + e _xe 12 ) + ... 

Central expansion 

1/2 

0(x,e) ~ h Q (x) + e h^fx) + eh^Cx) + ... 

Outer boundary layer expansion 

<Kx,e) - h 0 (l) - \ A 2 h 0 (l)(l-x) + e 1/2 h 1 (l) 

+ A 2 (A 2 -2)h 0 (l)(l-x) 2 /16 

- \ e 1/2 (A 2 h 1 (l) + A 2 h Q (l) ) (1-x) 

- (A 2 -6)(A 2 -2)A 2 h 0 (l)(l-x) 3 /288 
+ h 2 (l) + ... 

2 2 

In these expressions Cq» h^(x), h^(x), A^ and A^ come respectively from 
Equations (48), (21) and (42), (52), (43), and (50). Note that despite the 
large number of terms evaluated in the outer boundary layer the final expansion 
contains no terms involving powers of e above 1/2. 
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Diametral Beam Outer Expansions 


The differential equation and central solution expansion are the same 
for radial and diametral beams, so that with the exception of the boundary 
conditions in Equations (2) - (3) and (15) - (16) the development is the 
same for both beams through Equation (27), with k = 1 for the radial beam 
and k > 1 for the diametral beam. Different boundary conditions imply differ- 
ent boundary layer solutions, and in the case of the diametral beam symmetry 
indicates that the same boundary layer expansions apply at x = k""* - and x = 

In what follows we develop the solution near x « k \ 

As previously, we define a stretched coordinate 


X = <k 1 - xJB" 1 ^) 

where B(e) + 0 as e + 0, and Introduce the expansion 
<{)(x,e) = YqCOFqOO + YjCO^Cx) + ... 

where y^(e) constitute an asymptotic sequence, i = 0, 1, 2, 
(13) then becomes 


rY 0 (£) d 4 F 0 Y x (e) d\ 1 

L 3 4 (e) dx 4 B 4 (e)d^ + * ‘ * J 


- 1 - (k _1 - B(e)x) 2 j£ 


^2-ir Y o (e) d F 0 + d F i + 


B 2 (e) dx 2 


o —2 

6 (E) dx Z 


-i _ ry e > dF o . dF i 


- 2(k - B(e)x) 


3(c) dx 3(c) dx 


ip n 

!l + 

!x 


(94) 


Equation 


] 


- An [Yn< e ^ F n + Yi(e)F. + ...] 


- (e)A^[YQ(e)FQ + Y 1 (e)F 1 +...]- 0 


(95) 


Under the assumption that pretension T is sufficient to maintain 1 - k 2 = 0(1), 
we can choose 3 • 


3(0 - e 1/2 


(96) 
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and obtain from Equation (95) the dominant boundary layer equation 

,.*2 „ 

■™ K. a = U 




dx 


( 97 ) 


a 2 ““2 

where k A. (1 - k ). Equation (97) can be solved as was Equation (31), since 


Fq and gg have the same boundary conditions (see Equation (33)) with the result 

Fq(x) = A Q (k*x - 1 + e -k * x ) (98) 

where Aq is to be determined by matching with the central solution expansion. 
Matching requires the introduction of the intermediate limit defined by 

(99) 


x„ = (k 1 - x)n ^(e) 


1/2 

with r|(£) -*• 0 as e -*• 0 and e = o(n(£))» so that in the limit as e 0 with 


x^ fixed, (k 1 - x) = x^ri ->0 and x = (k - * - x)e 


- 1/2 - 1/2 

= x ne 7 ->-05, Matching 


in this limit requires, from Equations (18) and (95) 


Lim 

e+0 

x^ fixed 

n 


h Q (x) + v 1 (e)h 1 (x) + ... - Yq(e)F 0 (x) - Y 1 (e)F 1 (x) - 


-1 


Expanding in Taylor series about x = k and substituting Equation (98) 
produces 


Lim 

e+0 

x fixed 

n 


h 0 (k h - hg(k" 1 )nx ri +|hg(k 1 )n 2 xJ + ... + v 1 (e)h 1 (k" 1 ) 

- 1/2 

1 1 1 t ry — k^E tlX 

- v 1 (e)h 1 (k Jnx^ + ... - Yq(e)Aq ( k e ' nx - 1 + e ) 


- 1/2 

- Y 1 (s)F 1 (e nx^) - ... 


- o 


Matching requirements suggest the following: 
h 0 (k _1 ) = 0 
Y 0 (e) = e 1/2 

* 1 —I 

A Q k = -h Q (k X ) 


( 100 ) 

( 101 ) 

( 102 ) 

(103) 
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(105) 


A q = - h^k" 1 ) 

With this identification, we can proceed to establish h^(x) and F^(x), and 

then return to Equation (100) to confirm the matching. We can also use 

Equation (101) in conjunction with the Legendre function solution for h^ 

in Equation (21) to obtain values for A^. 

The singularity of the Legendre functions for unity argument does not 

present a problem for the diametral beam, since k" 1 < 1. Thus both U (x) and 

n 

V Q (x) are permissible solutions, subject to Equation (101). The symmetric 

or "even" modes are represented by U n (x), and the antisymmetric or "odd" modes 

are represented by V r (x), as indicated by Equations (22) and (23), which 

2 

together with Equation (101) provide solutions for A^. These results are 
illustrated numerically in the section which follows the present section on 
the boundary layer expansion. 

Orthogonality of the modal functions <f>(x) corresponding to different 
2 2 2 

discrete values of X , say X and X , can be established as for the radial 

(X OC p 

beam, to obtain the orthogonality relationship (as in Equation (45)) 
k -l 

4-1 * 5 cB < 106 > 

For the terms in the asymptotic expansion of a single modal function, we 
now have (in parallel with Equations (46) and (47)) 


-1 


f i = 1 

<k 


(107) 


k -l 

J\. i h Q (x)h 1 (x)dx = 0 


( 108 ) 


under the same conditions affecting the boundary layer solution. 



Thus we have for the even inodes, from Equations (107) and (21), with B = 0, 

/*k 2 2 

1-1 AV(x)dK = 1 




so that 


[ & < 


(x)dx 


-1 


(109) 


and similarly for the odd modes A = 0 and 

-1 „ -1 


B 




(x)dx 


( 110 ) 


Equations (107) and (108) indicate that Equations (49) and (27a) apply 

2 

to the diametral beam also. We can now obtain a new expression for by 

— 1 —1 

multiplying Equation (27a) by h Q and then integrating from -k to k , to 
obtain, after noting Equations (107) and (108), 

-1 


f k 

4 1 


f (l-x 2 )h^]’h 0 dx = -A 2 


Two integrations by parts, utilizing boundary conditions established by 
Equation (101) and its counterpart at x ■ -k \ and involving Equation (20) 
and Equation (108), produce 


A 


2 

1 


(1 - x 2 )h Q h 1 



For the odd modes the slope is the same at both ends but the displacement is 
opposite in sign; for the even modes this relationship is reversed. Thus in 
any case 


A 2 = 2(1 - k 2 )h 0 (k" 1 )h 1 (k~ 1 ) 


dll) 
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By using Equation (105) and then Equation (103), we can rewrite Equation 
(111) as 


.2 0,-1/, 2 n l/2 r 

A x = 2k (k - 1) [h 0 (k. )] 


(112) 


To obtain numerical answers from Equation (112) one must use the derivative 

of h^(x) as found in Equation (21) together with Equation (109) for even modes 

and Equation (110) for odd modes. 

2 

Now we could return with to Equation (27a) to solve for h^(x), using 
the boundary conditions established by matching in Equation (100) . Satisfied 
that h^(x) is fully determined, we now seek to establish the second term in 
the asymptotic expansion in Equation (95) . The matching requirement in 
Equation (100) suggests that 

Y 1 (e) “ e (113) 

Under this assumption, and with Equations (102) and (96), we can obtain from 
Equation (95) the following equation to be solved for F^(x): 


A 1 „ . -2, A 1 „-l dF 0 , A 0 

— - (1 - k > — * 2k w + 2k x tt 


dx 


dx 


dx 


or, in view of Equations (98) and the definition of k , 


1 ~ k * 2 = 2k _1 A 0 k*(l - e -k * X + xk*e k * x ] 


dx 


dx 


(114) 


Recall that A^k is fully determined by Equation (103) . 

Equation (114) has the same homogeneous solution as Equation (97), and in 
addition the particular solution F-^x) such that 


d 2 F 


dx 2 


1 r 1.1 k*X 1 —2 -k*x, 

2k A 0 k [ - - 5 fj + — xe - T x e ] 

k 4k 


(115) 


Proceeding again in parallel with the solution of Equation (31), we observe 
that for the total solution 
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A . -k^. d2F l 

— - A x e + -=2~ 
dx dx 

After two integrations, with integration constants A ^ and A^ introduced, and 
with A^ and A^ removed by imposition of boundary conditions on F^(x) and 
dF^/dx, the final solution becomes 


2k 


k- X A 


^x 2 + 


(__0 _ k ft A 3 j x + A. 


(116) 


Since Aq in Equation (116) is known from Equation (103) , the only free 
constant in this expression is A^; this unknown will be established by the 
matching process. Examination of Equation (116) reveals that matching in 
Equation (100) requires that 


Lira 

£+0 

x fixed 

n 


\ h^ (k“ 1 )n 2 x 2 n - v^Oh^k' 1 )nx n - Y 1 (e)F 1 (£' 1/2 nx ) 


= 0 


(117) 


If the first three terms in F 1 (x) in Equation (116) are ignored as trans- 
cendentally small, and Equations (104) and (113) are substituted, Equation 
(117) is satisfied up to terms below e if 


2k~ 1 A, 


*2 


0 * ’ -1 
“ - k A 3 = h x (k X ) 


and 


k 1 A 


TT 9 - “ " I ^Ck’ 1 ) 


(118) 


(119) 


Equation (118) provides the required value for A 3 , but Equation (119) 
must be an identity if the proposed solution is to be verified. This equation 
can be confirmed by using Equation (103) to rewrite Equation (119), noting the 
definition of k , as 
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( 120 ) 


(1 - k' 2 )hQ(k _l ) - 2k" l hg(k X ) = 0 

The validity of Equation (120) follows from Equation (20) with x •» k” 1 , since 
V k 1 > = 0 Equation (101). Thus the proposed solution is established. 
Diametral Beam Solution Summary 

The rotating pretensioned diametral beam in Figure 2, as characterized 

by transverse vibration equations (Equation (10)) with fixed end boundary 

conditions (Equation (16)) has natural frequencies given by Equation (17), 

2 

with X from Equation (19) and (49) available as 


,2 .2 1 / 2.2 . 

A - Aq + e Ai + 


2 

where ^ is given by Equation (112) in terms of h Q from Equation (21), and 
2 

where A Q is given for even modes by the combination of Equations (22) and (101) 
in the form 


u n (k ) 1 21 y 0 + 4! M 0 (fJ 0 6k 2 ' ) 

‘ y 0 Cy 0 " 6k " 2 >^0 “ 20k_2) + = 0 (12D 

2 2 

where Pq A A^/k , and for the odd modes from Equations (23) and (101) by 

V n (k_1) “ 1 - IT (U 0 " 2k ~ 2) + JT <4 - 2k ' 2 >^0 • l2k_2) - ••• - 0 

( 122 ) 

These equations have obvious solutions for particular values of k ; for 

example with p 2 = 6k -2 Equation (121) yields p 2 « 2 and hence k 2 = 3 and 
2 2-2 

Aq ° and with Pq = 20k Equation (121) yields either the '"first mode" with 

2 2 2 ? 

= 14.8 and hence k =1.35 and A Q = 20 or the "third" mode with Uq - 2.32 

2 2 

and hence k =8.6 and A^ = 20. More generally however one finds for any 

value of parameter k a spectrum of discrete solutions for p 2 . In the limiting 
2 -1 

case k + °° the Legendre function U n (k ) for the even modes becomes the 
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2 2 

cosine function, implying y^ “ tt /4, and corresponding to the vibrating taut 

2 9 

string, and in the neighborhood of the singularity at k * 1 we find y* ■+ 0. 

-1 2 

Similarly for the odd modes V^(k ) becomes the sine for k + ®, implying 
2 2 2 2 

Mq = V > and at k = 1, y^ = 2. Figures 4 and 5 portray the variations of y^ 

2 

(and hence natural frequencies) with k for the first and second modes 
respectively. 

In addition to the natural frequency equations, we have developed 
expressions for mode shapes in the form of asymptotic expansions as follows: 
Central expansion 

<Kx,e) * hg(x) + e 1 // 2 h 1 (x) + ... 

Boundary layer expansion near x = k ^ 


<Kx,e) 


1/2 ' -1. 

_1 t ^ h/} (k ) 

-Ck 1 - x)h 0 (k L ) + — — g n 


. e -k(k~ 1 - x ) E " ;L/2 


] 


-1 , - 1/2 


+ k 1 h’(k' 1 )k*(k 1 - x) 2 e k * (k x)e 


' „ -1 


+ 1 5«£J. a- 1 - ^(k-W 172 


k 


+ e 


k 1 


k^hjck" 1 ) 2 , , 

*2 Ck - x) + e ' h x (k )(k 1 - x) 


'2kV(k X ) , 1 

- £ [ -*4 + h 1 (k“- L ) + ••• 


and a similar boundary layer near x = -k” 1 . In these expressions k* 2 = (1 - k 2 ) , 
^(x) is the solution of Equation (27a) with A 2 from Equation (112), and h 0 (x) 
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is given by Equation (21), with the constants A and B in Equation (21) 
established by Equations (109) and (110) . Figures 6 and 7 are plots of 
h Q (s) for modes one and two respectively, with s A £/L. 

Conclusions 

The method of matched asymptotic expansions has a unique advantage in the 
modal analysis of rapidly rotating beams, in that the results are literal 
rather than numerical, and thus are applicable at once to a range of values 
of the system parameters, which include beam density, stiffness, pretension, 
and spin rate. The primary disadvantage is the approximate nature of the 
results, although the order of the approximation is well established. 

Application of this method to modal analysis of rotating plates is also 

12 

feasible in some cases, as will be reported in a separate paper. 
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